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Abstract 


The Bulk Lunar Electrical Conductivity 

by 

Donald Lucien Leavy 


Submitted to the Department of Earth and 
Planetary Sciences on January 29, 1975 
in partial fulfillment of the 
requirements for the degree of 
Doctor of Philosophy 


We study the electrical conductivity structure 
of a spherically layered moon consistent with the very 
low frequency magnetic data (0.0002 ^ f - 0.04 Hz) col- 
lected on the lunar surface and by Explorer 35. In order 
to obtain good agreement with the lunar surface magneto- 
meter observations, the inclusion of a void cavity behind 
the moon requires a conductivity at shallow depths higher 
than that of models having the solar wind impinging on 
all sides. By varying only the source parameters, a con- 
ductivity model can be found that yields a good fit to both 
the tangential response upstream and the radial response 
downstream. This model also satisfies the dark side tan- 
gential response in the frequency range above 0.006 Hz 
but the few data points presently available below this 
range do not seem to agree with the theory. 

A common feature of models resulting from the in- 
version of the sunlit side data is that the electrical 
conductivity profiles hardly increase by one order of 
magnitude at depths between about 200 and 700 km. Two 
simple interpretations of this constraint appear mutually 
exclusive at this point. On one hand, the 'persistence 
of a large temperature gradient to moderate depths, in 



models resulting from conventional thermal history cal- 
culations, would seem to require a conduction mechanism 
characterized by a very low activation energy {0.09 - 
0.24 ev) and moderate conductivity prefactor { 10“ 3 -10 -2 mho/m) 
On the other hand, the conductivity- temperature rela- 
tionships usually found in silicate minerals would lead 
to models of temperature -with very small gradients .at 
depths greater than about 200 km. 


Thesis Supervisor: T. R. Madden 

Title: Professor of Geophysics 
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Chapter I 

Historical Introduction 

1 . 1 Introduction 

Two major experiments have been the turning point, 
in recent years, in our understanding of the moon's 
electromagnetic environment. On July 22, 1967, the 
Explorer 35 spacecraft was injected into a stable lunar 
orbit carrying magnetometers, energetic particle detec- 
tors, and plasma probes on board. The results of this 
experiment were to unravel the essential feature of the 
interaction of the solar wind with the moon. However, 
it turned out that no effect of the conductive lunar in- 
terior could be detected at the satellite orbit. Such 
signals were not conclusively obtained until the deployment 
of the Apollo 12 lunar surface magnetometer (LSM) on 
November 19, 1969. Since the evolution of our concept has 
been largely shaped by the data obtained in these two ex- 
periments, they provide a natural division in the short 
history of the subject. 

1. 2 The Pre-Explorer 35 Period 

In order to develop a quantitative theory of the 
interaction of the solar wind with a planet, it is essential 
to know the magnitude of the steady, global magnetic field 
this planet might possess. At the earth's orbit, this 
field needed only to be of the order of 50 gammas in 
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order to balance the dynamic pressure of the solar wind 
and thus form a bow shock on the sunlit side of the 
moon (Willis, 1971). 

The results of the early spacecrafts sent to the 
moon were not entirely conclusive on this question. The 
Luna 2 probe (Dolginov, . 1961) did not observe any per- 
turbation of the interplanetary field at 50 km above the 
sunlit lunar surface, when the moon v;as in the magneto- 
sheath. The accuracy of the instrument was about 100 y, 
so this experiment did put a fairly accurate upper limit 
on the global magnetic field that might be present on the 
lunar surface. However, the possibility that a shock 
existed was not completely ruled out. A steady dipolar 
field of about 50 y presumably would be compressed to 
within a few plasma skin depths (=; 2 km) of the lunar 
surface and thus be hardly observable at an altitude of 
50 km . 

The interaction of the solar wind with an electrically 
conductive lunar interior might also build up the required 
50 y for a shock. Two modes of interaction are possible. 

In the poloidal H mode, eddy currents are generated in- 
side the moon by the time-varying interplanetary magnetic 
field. If we consider an homogeneous lunar model, these 
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currents confine the interplanetary magnetic field 

■ 

fluctuation to within a skin depth, 




(_? — ) 


1/2 


of the lunar surface .• Moreover, since the solar wind 
presumably shields itself from the magnetic field 
generated by these currents, through confining current 
flowing within a few plasma skin depths from the lunar 
surface, an interplanetary fluctuation associated with a 
lunar skin depth much smaller than the radius of the moon 
would be amplified roughly by a factor on the lunar 

surface. If we assume that a shock is produced by the 
sector .structure fluctuation of the 


interplanetary magnetic field, B ^ 5y , period s .10 days, 
(Schatten, 1971) , then we require a near surface lunar 
conductivity higher than a few mhos/m (6 0.1 R^) . 

Such high conductivity , typical of sea water on 
earth, was considered by Gold (1966) in a qualitative 
analysis of the interaction of the solar wind with the 


moon. Tozer et'al. (1967) quickly pointed out, however, 
that if water is not present inside the moon, the bulk 
lunar conductivity is likely to be determined by the tem- 
perature and composition at a given depth. Using a 
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temperature model proposed by Urey (1962) and the con- 
ductivity-temperature relationship for an olivine with 
10% fayalite, they showed that a conductivity of a few 
mhos/M is likely to be reached only deep inside the 
moon (R - 800 km) . They thus dismissed the possibility 
for the poloidal H mode to produce a detached bow shock 
in front of the moon. However, apparently omitting the 
possibility for the magnetic field line to slip around 
the core, Tozer et al . concluded that the field lines will accrete 
in front of the core and thus produce an attached hydromagnetic 
shock near the limb of the optical shadow. 

The toroidal H interaction was investigated by 
Sonett et al . (1967). In the rest frame of the moon, 

the interplanetary magnetic field, Bg^, is seen to be ac- 
companied by an electric field given by 





x B 


SW 


( 1 . 2 . 1 ) 


An exact determination of the solar wind velocity, 
V-..., involves not only the streaming speed of the 

oW 

solar wind (- 350 km/sec) , but also the various motions of 
the moon.. However, even the most important of these motions, 
the rotation around the sun together with the earth 
(- 29.8 km/sec) , has a magnitude much smaller than the 
bulk solar wind velocity. Thus these motions wall be neglected 
in the following discussion. 
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Let us approximate the moon by a cylinder with 
axis parallel to the E field, with radius equal to the 
moon radius and with length twice this value. The po- 
tential across the cylinder is then given by 


4> 2R m E 2 R m V sw B SW 


( 1 . 2 . 2 ) 


where we have assumed the solar wind magnetic field to 
be perpendicular to the solar wind velocity. 

Equation (1.2.2) also assumes that the electric con- 
ductivity of the solar wind is very high so that no 
electric field is seen in the rest frame of the wind. 
By Ohm * s Law , we have 


Re 2 ^ 


(1.2.3) 


where I is the current and Re the resistance along the 
cylinder axis and cr is the cylinder conductivity. From 
Ampere's Law, and combining Equations (1.2.2) and (1.2.3) 
we obtain 


B 


toroidal 


Hd 1 

2 ttR m ay o aR M V SW B SW 


(1.2.4) 


where for a homogeneous cylinder a = 0.5. 

We note that, for field fluctuations associated with 
an interplanetary wavelength and lunar skin depth much 
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larger than the radius of the moon, the toroidal H in- 
duction field is independent of frequency. Moreover, if 
. . . -5 

the conductivity is about 3 x 10 nho/m , the toroidal H 
field at the lunar surface is nearly ten times the solar 
wind magnetic field. In that case, a shock might be 
formed on the sunlit side of the moon. However, even if 
the interior conductivity of the moon is higher than 
10 mho/m but is covered by an insulating layer, the 
toroidal H field may become much smaller than the solar 
wind field. Let us assume our cylinder to be capped by 
layers of thicknesses A and neglect the internal re- 
sistance compared to the one at the surface (the resistance 
to the current flow added in series) . Then, the toroidal 
H field is still given approximately by Equation (1.2.4) 

but with a = 0.5 r^/A and a equal to the crustal conduc- 

-7 

tivity. A 17 km crust with conductivity 3 x 10 mho/m 
still gives a surface field ten times higher than the 
solar wind field. However, if the conductivity-thickness 

ratio of the surface layer is a hundred times smaller 

— 8 

than 2 x 10 mho/m-km, then the toroidal H field becomes 
only one tenth of the solar wind field-, at the surface of 
--the moon. 

In situ , measurement of rock conductivity at the sur- 
face of the earth has not revealed conductivity much 


OF POOR qPALH* 
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-5 *- 

lower than 10 mho/m. But this is due mainly to the 
presence of water in the earth crust (Madden, 1971, and 

Brace, 1971). Laboratory measurement on dry rocks have 

— 8 

shown, however, that conductivity lower than 10 mho/m 
can easily be reached at room temperature (see, for ex- 
ample, Fensler, 1962). Because of the relatively large 
range of conductivity one might assume for the surface 
of a planet, predicting the importance of the toroidal 
H mode will probably remain difficult. 

During a month in 1966, the Luna 10 spacecraft 
was placed in lunar orbit and sent back to earth additional 
data on the magnetic field around the moon. A very regular 
field of 23 to 40 y was observed (Dolginov et al . 1966) . 
This field did not vary much either along an orbit of the 
satellite around the moon (periselene: 1.2 R^; aposelene: 
1.7 R^) or along the orbit of the moon around the earth. 
Though the regular behavior of the field rendered the 
measurement somewhat suspect, Dolginov et al . (1967) pro- 

posed a possible interpretation in terms of the inter- 
action of the solar wind^ with a conductive lunar interior. 

1 . 3 From Explorer 35 to Apollo 12 

Explorer 35 has a stable orbit of period 11.5 hours, 
aposelene = 5.4 R^, and periselene - 1.4 R^. The pre- 
liminary results sent back to earth by this spacecraft 
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(Ness et al . , 1967 , and Lyon et al . , 1967) were to be 
confirmed by subsequent instruments sent to the moon. 

Contrary to what was observed by Luna 10, no 
steady lunar magnetic field, of magnitude several times 
the average solar wind field, was found. Behannon (1968) 
by examining the Explorer 35 data, obtained when the moon 
traversed the neutral sheet of the geomagnetic tail, was 
able to distinguish the magnetic field of a possible per- 
manent lunar dipole from the one induced by a bulk lunar 
permeability. He was thus able to establish an upper limit 
of 4 y on the permanent dipole field at the lunar surface 
and an upper limit of 1.8 for the bulk relative magnetic 
permeability of the moon. 

No bow shock wave was observed. As we have seen 
above, a possible explanation for the absence of a toroid- 
al H mode-induced shock is the presence of a surface layer 

of conductivity thickness ratio smaller than 

— 8 

2 x 10 mho/m - km. The fact that the time variation of 
the interplanetary magnetic field associated with its 
sector structure does not produce a poloidal H type of 
shock tends to imply that the conductivity in the top 
200 km of the moon is smaller than a fev; mhos/m. 

When the moon is in the solar wind and the satel- 
lite passes through a cylinder approximately defined by 
the optical shadow, several effects of the interaction of 


Original page ts 
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the solar wind with the moon were observed. The plasma 
flux (50 < Ep < 2850 ev) decreased by several orders of 
magnitude consistent with the hypothesis that the particles . 

in the solar wind were absorbed on the sunlit surface of 

* 

the moon, leaving a void on its downstream side. The 
magnetic field was also perturbed (Figure 1.1). Near the 
boundary of the optical shadow a small decrease in the field 
magnitude is followed by a gradual increase in the plasma 
umbra (Colburn et al . , 1967, and Taylor et al . , 1968). 
Several theoretical attempts were made to explain these 
characteristics. Spreiter et al. , (1970) , using the 

equations of magnetohydrodynamics, examined the case where 
the interplanetary magnetic field is aligned with the 
flow direction. This particular field configuration was 
reported by Ness et al . (1968) , and their results agree 

qualitatively with the observations. They found that the 
small decrease can be understood in terms of the approximate 
conservation of the ratio of magnetic field to particle 
density along a streamline. The solar wind particles 
initially tend to fill the void thus decreasing particle 
density, pressure and magnetic field along a streamline. 
These gradients are accompanied by current that increases 
the field of the plasma umbra. An equilibrium void/plasma 
boundary is reached when the umbral magnetic pressure 
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balances the particle and magnetic pressure of the solar 
wind. Whang (19 70) , noting that the scale length as- 
sociated with the observed field gradient is generally 
larger than the proton gyroradius 100 km) , used the 
guiding center approximation to calculate the field and 
particle distribution. His treatment included the effect 
of anisotropic propagation of the magnetosonic wave in 
addition to the more general field configuration. 




EXPLORER 35 


5 AUGUST 1967 


Figure 1.1 

Simultaneous measurements of field and 
plasma obtained on August 5, 1967, from lunar orbit 
on the Explorer-35 spacecraft. The trajectory of the 
spacecraft is shown projected on the ecliptic plane and 
positionally correlated with the data through UT annota- 
tion. The x axis is parallel to the sun-moon line. 
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A fundamental aspect of the interaction is the fact 
that the Alfven and sound speed 30 km/sec) are much 
lower than. the bulk solar wind velocity. Consequently, 
a given position in the solar wind, on the downstream side 
of the moon, will be influenced by the* disturbed con- 
dition at the plasma -vacuum boundary only if it is 
within the Mach forecones originating at the terminator 
(the limb for solar wind flow) . The existence of such 
Mach cone received some experimental support by Whang 
et al . (1970). Also, Ogilvie et al . (1970) showed that 
the amplitude of the umbral increase and penumbral de- 
crease tend to grow in proportion to the ratio of par- 
ticle to magnetic pressure. However, since, to date, only 

i-' ‘ 

a subset of the parameter that influences the field charac 
teristics have been compared to the theory, its detail- 
ed confirmation must still be considered incomplete. 

In addition to the two characteristics just men- 
tioned, a small increase of the magnetic field was often 
observed on the interplanetary side of the penumbral 
decrease. This feature was shown to correlate with a 
small enhancement in particle density (Figure 1.1). 

(See Siscoe, 1969). Hollweg (1968) suggested that if 
ice is present near the surface of the moon a significant 
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toroidal H field, might be induced. The interaction of the 
solar wind with this field would presumably produce the 
small external increase. However, chemical analysis of 
lunar samples brought back to earth in the Apollo 11 mis- 
sion did not show any evidence of hydrous phases in lunar 
rocks (see, for example, Charles et al. , 1971). Schwartz 
et al. , (1969) , pointed out, though, that a dry moon, with 

a hot interior, might also produce a significant toroidal 
H field, depending on its near-surface conductivity. But 
Ness (1972) , in a review of the subject, argued against 
this possibility on the basis that the small peak is often 
observed when the interplanetary magnetic field is aligned 
with the solar wind flow velocity. According to 
Equation (1.2.1), no motional electric field should exist 
in that case. 

In view of the absence of bow shock and in 
anticipation of the lunar surface magnetometer experiment, 
the response of the moon to time varying magnetic fluc- 
tuation was also re-evaluated theoretically by several 
workers (Blank et al . , 1969, Schubert et al. , 1969, 

Schwartz et al . , 1960, Sill et al . , 1970). They showed 
that the expected frequency dependence of the poloidal H 
lunar response might enable us to distinguish between 
several possible lunar conductivity profiles, in particular 
between a hot and cold moon based on the conductivity 
models proposed by England et al. (1968). 
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The concept that the solar wind interaction with the 
moon does not perturb regions outside the Mach cones on 
the downstream side of the moon was found to break down • 
for frequencies higher than approximately 0.1 Hz (Ness 
et al . , 1969) . High frequency fluctuations seemed to 
originate at the plasma vacuum interface and to propagate, 
outside the plasma umbra, along the time average magnetic 
field lines threading the lunar wake. A higher level in 
the power spectral density above about 0.1 Hz was seen 
to be maintained in a region within approximately one 
lunar radius of the plasma/void interface. A tentative 
explanation for this phenomenon was given by Krall et al . 
(1968) in terms of an electron ballistic effect. They 
argue that the solar wind electrons (thermal-speed 
= 2000 km/sec, gyroradius - 2 km) , upon reflection at the 
solar wind/void interface, will carry a memory of the 
perturbed condition at this boundary. This memory, which 
manifests itself as a fluctuating magnetic field as- 
sociated with current produced by a perturbation of the 
electron distribution function, will eventually fade away 
by phase mixing as the electrons travel away from the 
boundary. 

1 . 4 From Apollo 12 to Recent Years 

Several magnetometers were flown to the moon during the 
Apollo missions. One highlight of these experiments was 
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the discovery of local remanent magnetic fields at most 
Apollo sites (see Table 1.1 from Dyal et al . , 1972, 1973 a 
b). 


Site 

Apollo 12 
Apollo 14 

Apollo 15 
Apollo 16 


Table 1.1 

Location Steady Magnetic 

Field (y) 


Oceanus Procellarum 38±3 

( 3 . 2°S 23 . 4°W) 


Fra Mauro 
(3 . 7°S 17 . 5°S) 

(two sites 1.1 km 
apart) 

Hadly-Apennines 
(26.1°N 3.7°E) 

Descartes 
( 8 . 90s 15 . 5°E) 

(five sites separated 
by 0.5 to 7.1 km) 


103 + 5 
43±6 

6±4 

327 
232 
1.8 9 
• 113 
113 




In addition to these instruments, two subsatellites, 

with magnetometers on board, were launched during the 

Apollo 15 and 16 missions along orbit near the surface 

of the moon (^ 100 krn).. Data from these subsatellites 

were used to map the lunar remanent field (Coleman et al . , 

18 

1972) and also to place an upper limit of 3.6 x 10 gauss 
3 

cm on the permanent magnetic dipole moment of the moon 


(Russell et al. , 1973). Thus, the permanent lunar mag- 
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“8 

netic moment is at most 4.5 x 10 times as strong as the 
one of the earth and can only produce a surface field 
smaller than 0.2 y. A review of the magnetic proper- 
ties of lunar rocks can be found in Fuller (1974) . 

New bounds on the bulk relative magnetic permeabi- 
lity of the moon (— = 1.03±0.02) were established by 

^ o 

Parkin et al . (1973) by considering the data from 

Explorer 35 and the lunar surface magnetometers when the 
moon was in the geomagnetic tail. 

No significant toroidal H field was found on the 
surface of the moon. However, the discovery of a large 
remanent magnetic field leads to the suggestion by Barnes 
et al . (1971) that when such a field is present near the 

terminator, it might interact with the solar wind to 
produce the observed small limb compression. A prelim- 
inary check of this hypothesis was made by Lichtenstein 
et al. (1974) with the data from the Apollo 15 subsatel- 
lite. They found that the occurrence rates of limb com- 
pression are roughly proportional to the amplitudes of 
the remanent fields observed at satellite altitude. 

Large amplification of the tangential components of 
the magnetic field at the lunar surface was observed by 
comparing the magnetic fluctuations at Explorer 35 to 
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the one at a LSM when both instruments were on the sunlit 
side of the moon and outside the geomagnetic tail (see, 
for example, Sonett et al . , 1971 a). Detailed power spectral 
density of each component of the field was evaluated at the 

T T 

LSM (P TCT ..) and Explorer 35 (P„ vr> } when the latter instru- 
ment was bn the upstream side of the moon. If, at an LSM 
site, we define a mutually orthogonal set of unit vectors 

_y ■ y 

x, y, z (vertical, eastward and northward, respectively), 
then the data can be expressed in the form of an ampli- 
fication factor A. (L - x, y, z) were 



(P L /P L ) 
v LSM' EXP y 
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The coordinate system used to evaluate the power 
spectral density from Explorer 35 is made to correspond 
to the one used at the LSM site. The data- is presented 
in Figures 1.2, 1.3, and 1.4. These power spectra were 
evaluated when the moon was either in the free Streaming 
solar wind or in the earth's magnetosheath. Figure 1.2 
represents data obtained at the Apollo 12 site when the 
LSM was on the sunlit side of the moon. The data in 
Figure 1.3 were also obtained at the Apollo 12 site but the 
LSM was on the night side of the moon, within 45° from 
the antisolar point. Figure 1.4 represents data ob- 
tained at the Apollo 15 site when the LSM was on the 
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sunlit side of the moon, within 45° from the terminator. 
The data in Figures 1.2 and 1.3 are from Smith et al. 
(1973), whereas those in Figure 1.4 are from Schubert 
et al. (1973) . The technique used to evaluate these 
spectra has been discussed by Sonett et al . (1971 b) . 

Recently, Schubert et al. (1974 b) have published the 
amplification factors obtained when the moon was in the 
plasma sheet of the geomagnetic tail- These latter data 
will not be used in this thesis. 

In the initial inversion of the sunlit side data, 
Sonett et al . (1971 c) assume the solar wind plasma to 

be radially incident on the moon. They obtain a conduc- 
tivity profile characterized by a relatively large spike 
at a depth of approximately 250 km. But conductivity pro- 
files with this peaked behavior were soon recognized to be 
only part of a larger set of models, many of which 
smoothly varying, that fit equally well the frontside 
data. (See, for example, Kuckes, 1971 and Sonett, 1972.) 

The night side data were initially interpreted by 
Dval et al . (1973 b }, who investigated the passage of 

large magnetic field discontinuities in the solar wind. 
They assumed the moon to be surrounded by a vacuum and 
inverted the time domain response of the moon to these 



discontinuities.. They found that, in addition to a 

constant conductivity layer of 3 x 10 -4 mho/m in the 

upper 700 km of the moon (except for a thin — - 40 km — 

non-conducting crust) , they also require a core of con- 
, , -2 

ductivity 10 mho/m at depth greater than 700 km in 
order to fit the decay time of the vertical magnetic field 
component of the discontinuity. Schubert et al. (1973 a) , 
however, showed that, using the vacuum approximation, the 
radial amplification factor in Figure 3 poorly,, resolves 
the conductivity of the bottom layer. Moreover, they 
showed that the tangential amplification factors on the 
night side cannot be inverted within a model that assumed 
the moon to be embedded in a vacuum. This is due to the 
neglect of confining current on the frontside of the moon 
and at the plasma vacuum boundary. The effect of these 
currents is to amplify the tangential surface magnetic 
field to a value higher than can possibly be reached by 
any conductivity model within the vacuum approximation. 

The radially incident plasma model also suffers from 
its neglect of the plasma void behind the moon. One 
aim of this study is to incorporate in a single model the 
dayside-nightside asymmetry in the plasma environment of 
the moon. Concurrent with our effort, Schubert et al . 


0.973, b, c) , Schwartz et al. (1973) and Smith et al. 
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(1973) have recently published initial results of a 
theory which account for the day-nightside asymmetry. 

We shall discuss some of their results as we go along 
in the next chapters. 

1 - 5 Thesis Outline 

Chapter II contains the field solution in the 
different regions around the moon, together with a dis- 
cussion of the boundary conditions used in the analysis. 

In Chapter III we discuss the numerical method used 
to solve the forward problem together with the properties 
of the solution for various lunar conductivity models and 
solar wind parameters. 

In Chapter IV we solve the inverse problem for 
particular sets of parameters of the source field and 
discuss the resolvable feature of the conductivity 
structure together with the information distribution 
among the observations. Also, we subject the moon to 
various temperature models and discuss inferences that can 
be made on the parameters of a semi-conductor satisfying 
both the thermal and magnetic constraints. 
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Chapter ±1 

Analysis of the Field 
2 . 1 Introduction 

We will consider only the period of the lunar month 
when the moon is either in the earth’s magnetosheath or 
in the free streaming solar wind. During this period 
we can distinguish three regions in the moon’s electro- 
magnetic environment: the solar wind, the void 

cavity behind the moon, and the interior of the moon. 

In the following sections we discuss the field repre- 
sentation in each of these regions and their coupling 
through the boundary conditions . 

The different coordinate systems used to represent 
the field are shown in Figure 2.1. Their origins should 
all coincide with the center of the moon but for clarity 
they have been translated parallel to the solar wind 
velocity. We will follow Morse et al. (1953) for the 
definition and notation of the various functions used 
in the text. 

Since we expect the lunar response to be somewhere 
between the one expected for a spherically symmetric plasma 
and a vacuum environment, it is instructive to compare 
some of their characteristics before attempting to in- 
corporate the various regions around the moon into a 
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more comprehensive model. This is especially true in 
our case since the solution of the more realistic model 
can be reached only through a numerical method which 
lacks, to some extent, the insight provided by a simple 
close-form formula. 

( 

To avoid any of the complications arising from the 
structure of the source field, we shall assume the inci- 
dent field to be spacially homogeneous and given (see 
Figure 2.1) by: 

H inc =a h a = H (a sin6sin0+a Q cos9 sint|)+a , cos4s ) (2.1.1) 

oy or 0 <p 

In the vacuum approximation, we allow the 
induced field to expand into a void outside the moon 
We can express the total field in the void as the sum of the in- 
cident field and the field of a dipole, i.e. : 

r m 

H vac = H lnC + V x V x b — sinBsincfa (2.1.2) 

r T r 

where r is the distance to the center of the moon and 
b is a constant to be determined by the boundary con- 
dition. 

In order to gain some insight into the amplitude 
of the response, let us assume that the moon is formed 
of two layers, the top one an insulator and the bottom 
one a perfect conductor. The field in the insulating 
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shell is given by 


gmoon ■ = y 


x V x 


4 2 . - 

(c — + dr ) sin0sin4>a r 


(2.1.3) 


In both models, the magnetic field inside the moon 
must have a vanishing normal component at the surface of 
the perfect conductor of radius a_. In the vacuum ap- 
proximation, the three components of the magnetic field 
must be continuous at the lunar surface. In the sym- 
metric plasma assumption, however, since we assume the 
induced field to be confined by a thin current layer 
above the surface of the moon, it is only necessary to 
equate, at the lunar surface, the normal component of the 
magnetic field inside the moon to the one of the incident 


field, as given in Equation (2.1.1). If we apply these 
boundary conditions and extract from the result the ratio 
of the tangential component of the magnetic field at the 
lunar surface to the one contained in the incident field, 
we obtain 


2 4 


3 a 3 

A plasma = 1 + 2 (2.1.5) 

R,, - a J 
M • 

These responses are plotted in Figure 2.2. We note 
first that the vacuum response can only reach an upper 
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Figure 2.2 
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limit of 1.5. Since the data of Figure 1.3 exhibits 
measured value higher than this when the LSM is near the 
antisolar point, one can readily infer that the vacuum 
approximation is inadequate to interpret the tangential 
amplification factor in the more realistic geometry. 
Another point of interest is the detectability of con- 
ductivity deep in the lunar interior. For example, if 
we assume the data to have a 10% relative error due to 
noise, we note that even a perfect conductor of radius 
lying between about 700 km (a - 0.4R^) and 1000 km 
(a = 0.6 R^) , would yield a response only of the order of the 
uncertainties in the data. 

A similar treatment can be given for the toroidal H 
magnetic field but since these results are completely 
worked out in Sill (1970) and are similar to our deriva- 
tion in a simpler geometry given in Section 1.2, we 
omit them here and proceed directly to the task of 
representing the field in the various regions of the 
lunar environment. 

2 . 2 Representation of the Field in the Solar Wind 

In addition to its quasi-stationary sector 
structure, the interplanetary magnetic field is usually 
permeated by various magnetohydrodynamic shocks, dis- 
continuities and waves. (For a review of the first 
two types of disturbance, see Burlaga, 1972). Belcher 
et al. (1971) have found substantial evidence that the 
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power spectral density of the interplanetary magnetic 
field, in the frequency range where lunar induction 
has been measured (Figure 1.2), is dominated at least , 

, fifty percent of the time by large amplitude Alfven 
waves, propagating outward from the sun. Sari et al. 
(1969) have proposed an alternative model to explain 
these micro-scale (< 0.01 au) fluctuations. They suggest 
that the interplanetary magnetic field has a spaghetti- 
like filamentary structure, with bundles of lines-of- 
force separated from one another by tangential discon- 
tinuities and essentially static in the rest frame of 
the solar wind. Though the physics involved in these 
two models differs markedly, these differences are of 
little consequence for our application. A common charac- 
teristic that must be retained, however, is that the inter- 
planetary magnetic fluctuations, as seen in the rest 
frame of the moon, are essentially convected at the solar 
wind speed. This approximation seems appropriate for the 
free streaming solar wind but the increase in the plasma 
pressure in the magnetosheath might produce wave velo - 
cities which are a substantial fraction of the solar 

r , 4 

wind velocity. Unfortunately, the data presented in 
Chapter I did not differentiate between these two lunar 
environments, so some caution should be used in their 
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interpretation. We' should note here that for 
waves propagating nearly perpendicular to the solar wind 
velocity, their phase velocity plays an important 
role in determining their spatial structure. Since 
such waves probably represent only a small fraction of 
the power density spectrum of the fluctuation, their con- 
tribution to the response shall not be considered any further. 

We shall assume that we can represent the fluctuation 
as a superposition of plane waves convected as a constant 
solar wind velocity, V . In the first of these hypotheses, we 
assume that the magnetic fluctuations are not generated 
near the moon and thus have little sphericity in their 
structures. In the second assumption we ignore perturbation 
such as large solar-flare-associated 
shocks which significantly modify the solar wind velocity 
A change of 100 km/sec in the bulk speed within a 
minute's interval can occur during such events (Chao, 

1970) . 

Let us 1 choose the y-axis in Figure 2.1 such that 
the wave normal to a given fluctuation is in the y-z 
plane and subtends an angle ip with respect to the solar wind velocity. 
If we assume ip to be positive if the wave normal has a 
component along the +y direction, we can express the 
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field as follows: 


H = H (a cosii-a sinij/)e 
o y r z y 


ik ( zcos^+ysimjj) 


+ H a e 
ox x 


ik (z cosi|;+ysin^) 


where 


kcosip - k // 


0 ) 


SW 


( 2 . 2 . 1 ) 


ksinty = k | = 


0 ) 


1~ v 


tanij; 


SW 


and 

* - -»o 5 sw * » 

A time factor, e is implicitly used here, 

as in all similar, equations of this chapter, and w 
stands for the angular frequency measured in the rest 
frame of the moon. In the solar wind rest frame, however, 
a fluctuation with phase speed V^ h and propagating paral- 
lel to the solar wind velocity is seen with an angular 

ojV h 

frequency equal to — — — . 

V SW ph 

For example, an Alfven wave, at the highest angular 
frequency for which the lunar response has been measured, 
0.25 rad. /sec. (Figure 1.2), has an angular frequency of 
0.02 - 0.05 rad. /sec. in the solar wind frame. This is still 
about an order of magnitude smaller than the nominal 



proton gyrof requency in the solar wind (0.5 rad./sec) . Also 
for the case of an Alfven wave, we have neglected the 
electric field in the rest frame of the wind. Equation 
(2.2.1) should still give a good first order approximation, 
however, since this electric field is smaller than the one associ- 
ated with the motion of the solar wind by the factor: 



(5 - 10 ) . 


2 . 3 The Field Inside the Moo n 

We assume the moon to consist of spherical layers 
of constant thicknesses and homogeneous conductivities. 
Such a structure is in harmony with a lunar model where 
the conductivity is mainly determined by the composition 
and temperature and where both can be assumed to vary 
only along the radius. If thermal conduction is the 
main heat transfer process inside the moon, this would 
seem to be a reasonable assumption for the temperature 
distribution. However, if solid state convection plays 
an important role in heat transfer, Turcotte et al . 

(1972) have suggested that a significant asymmetry in 
the temperature distribution could be introduced at 
depth where convection occurs. 

Also, if extensive lateral inhomogeneities in a 
global radially-varying moon structure are located near 
a given magnetometer, they might significantly influence 
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the signal observed at this instrument# especially at the 
higher frequencies. Even though power spectral data have 
been published from only the Apollo 12 and 15 site 
magnetometers, already there is substantial evidence 
that a regional signature is detected at the latter 
site. Schubert ot al . (19 74c) , by examining the power 

density spectrum at different angles in a plane tangent 
to the surface at the Apollo 15 site, found that the distribution 
of power is strongly peaked along the northwest-southeast line 
at frequencies above approximately 5 mHz. This 

1 + 

direction of polarization seems to be observed consistently 
and is independent of the directional character of 
the solar wind power spectrum and of the position of the 
magnetometer in the asymmetric plasma environment. A 
likely expalantion for the anisotropic character of the response 
is that' some regional - inhomogeneities influence the data at 
the site. . Maxima in the power spectrum of the tangential magnetic 
field components were also observed at the Apollo 12 site. But since 
between approximately 0.001 and 0.02 Hz, the peak tends to 
align along a line parallel to the tangential component 
of the renanent field at this site (= 64° south of east) , 
it was tentatively attributed to a modulation of this 
field by fluctuations in the dynamic pressure of the solar 
wind (Sonett et al. , 1972). We examine briefly 
this suggestion in Appendix I, but it should be pointed 



out here that anisotropy in the power spectrum of the 
tangential magnetic field component is also observed by 
the Apollo 12 magnetometer when this instrument is on 
the dark hemisphere of the moon. Its detailed pro- 
perties have not yet been presented in the literature, 
but its existence can be inferred from Figure 1.4. So 
a regional influence on the data might also accompany 
the noise due to the remanent field at this site. 

Thus, even if the bulk of the moon can be model- 
led approximately by a radially-varying structure, we never 
theless face the problem. of extracting its global in- 
duction signal from data contaminated by noise and 
regional inhomogeneities. A first empirical step in 
that direction was made by Sonett et al . (1972) who 

estimated the values of the tangential amplification 
factor at the Apollo 12 site, along a direction such 
that this quantity is minimal (e.g., between 0.001 
and G.02, along a line parallel to ~E25°N, i.e., ortho- 
gonal to the direction of the tangential remanent field) . 

We show these data, A m j_ n ' i n Figure 2.3, together with 
the data of Figure 1.4 (and a data point from Figure 1.2). 
We' note that the Apollo 15 data, though probably in- 
fluenced by the cause of the peak along the northwest- 
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southeast line, agree fairly well with these A « 

mm 

values. This set of data thus probably represents a 
good first approximation of the global response of the 
moon. But, of course, in order to qualify this in- 
ference, a more detailed data analysis is required,- in 
conjunction with theories that account for the possible 
regional influences and other sources of noise. 

We now turn to the solution of the 
Maxwell equation inside a layered sphere. These solutions 
are well known so we shall only summarize 
some of the main results here. 

The continuity of the tangential components of the 
3 and 3 fields* can be applied to each of the two electro- 
magnetic modes independently at an internal boundary of 
the sphere. This permits us to determine, for each har- 
monic of each mode, the ratio of electric to magnetic 
field at the lunar surface. Instead of dealing with im- 
pedances, however, we found it convenient to define, in 
Appendix II, some related quantities, and T^, which per 
mit us to write the field as follows: 

* 

Here we assume the moon to have a constant permeability 
equal to the free space value. 
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where k = i wu o , (a is the conductivity of the sur- 

s OS s 

face layer) and V is defined in Equation (2.2.1). 
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The coefficients a Q w and b a ^ are the coefficients 

nm nm 

of the poloidal and toroidal H mode, respectively, and 
are to be determined by applying the boundary conditions 
at the lunar surface. 

2 . 4 The Field in the Void Region on the Downstream Side 

of the Moon 

In order to obtain a tractable representation of 
the field in the void region on the dark side of the 
moon, some assumptions must be made on the parameters 
of the solar wind. We shall adopt the following three 
hypotheses: 

1. We can treat the solar wind as a cold plasma. 

In other words, we shall neglect any effect of the ther- 
mal pressure compared to the magnetic pressure. 

2. We assume that the ratio of the solar wind 
bulk speed to hydromagnetic wave speed can be consider- 
ed infinite. 

3. We assume the solar wind velocity to be time- 
and space-independent both in magnitude and direction. 

We hasten to point out, however , that the solar 
wind 6 is on average equal to about 1, and finite 
8 effects are readily observable in terms of penumbral 
decrease and umbral increase of the magnetic field (Ogil- 
vie al . # 1969). Moreover, a finite, but high, Mach 
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number effect has been observed by Whang et al . (1970) , 

who measured a Mach cone angle of 8°. But, unfortunately, 
it is difficult to relax these assumptions or to calculate 
their probable influence on the data, in the frequency 
range where lunar induction is important, since (i) simul- 
taneous data of all relevant solar wind parameters are 
generally not available, and (ii) we lack a proper theory 
that accounts for these effects. (To date, even steady 
state theory predicts some parameter dependence of the 
magnetic field that has no counterpart in the observation 
(see, for example, Ogilvie et al . , 1969 ) . 

We shall content ourselves in this thesis to check 
that a theory that incorporates these assumptions does 
in fact reproduce some of the major frequency-dependent 
characteristics of the field observed around the moon. 

This procedure is somewhat unsatisfactory, however, since 
the agreement comes largely as the result of a success- 
ful search for a conductivity profile that satisfies the 
observations within the framework of the approximation 
theory. But the very fact that such a model can be 
found that agrees reasonably well with both front and 
backside data, does give some measure of confidence in 
the approximation theory used. This trust is further 
improved by the fact that the theory can reproduce at 
this stage, at least qualitatively, some of the major 
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frequency characteristics of the field observed in the 
void region, sufficiently far downstream so that it is 
uninfluenced by the moon. 

In order to discuss this last point, let us ex- 
amine some of the consequences of the hypothesis we 

j ~ 

iiiauc . 

First, the geometry of the vacuum region can 
be modelled by a semi-infinite circular cylinder with its 
radius equal to the lunar radius. Moreover, the 
boundary conditions on the electromagnetic field can be 
most easily obtained since the boundary layer between 
the void and the undisturbed solar wind can be considered 
infinitesimally thin. This permits us to natch, at the 
boundary of the cylinder, the normal component of the 
solar wind magnetic field and its wavelength parallel 
to the axis of the cylinder to a solution of Maxwell’s 
equations inside the void. This solution can be given 
in terms of a cylindrical TE mode that also satisfies the 
continuity of the tangential electric field components 
across the plasma-vacuum boundary layer. The normal com- 
ponent of H of the solar wind field, at the boundary of 
the cylinder, can be evaluated from Equation (2.2.1) and 
the resulting field inside the void can be expressed as 


follows : 
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where 

e = Neuman factor = 1 when m = 0 and 
m 

= 2 when m > 0 

and where here and in the following development, we use 
the following definition : 
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i 

+ ... 

is odd 

= 0 

if 

m 

+ n 

+ 

i 

+ ... 

is even 


and 


0 if m + n + i + ... is odd 

1 if m + n + i + ... is even 

The prime stand for the derivative with respect 

to the argument and J , I are the Bessel and hyperbolic 
3 m m 

Bessel functions in their usual notation. 


Y 


mni . 


This field is not sufficient, however, to match 
the boundary conditions on the surface of the moon. To 
accomplish this, we must add an "end effect field" that 
will be discussed shortly. But before doing this, let 
us digress a little to examine some of the properties 
of Equation (2.4.1). 

We note first that for very low frequencies, 
ky/R^ << 1, the field is homogeneous and has the same 
value as the unperturbed solar wind field. Thus, due 
to our idealization of the plasma parameters, we do not 
reproduce the penumbral decrease nor umbral decrease 
observed behind the moon (Figure 1.1). On the other 
hand, for frequencies such that >> 1, Equation (2.4.1) 

shows that the pov/er spectral density in the void should 


be at a much lower level than what is observed in the 
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unperturbed solar wind. The ratio of the power spectral 
densities in that frequency range and on the axis of the' 

-2k, yK, 

cylinder, is approximately given by e /y/ . This charac- 

teristic of the field inside the void is due to the sur- 
face wave nature of the TE mode. The solution of the 
vector Laplace equation, propagating along the axis of 
the cylinder, decays approximately exponentially away 
from the boundary of the cylinder when kyR^>>l. In 
harmony with this result, a sharp drop in the power spec- 
tral density of the magnetic field in the void was observed 
by Ness e t al . (1969), at frequencies above about 0.1 Hz 

(k^yR^ - 3). But a detailed quantitative comparison is 
hindered by our lack of knowledge of the exact position 
of the satellite (in particular its distance from the 
axis of the cylinder) during these observations. 

Another type of observation that incorporates both 
of these limiting features is the response of the void 
cavity to large discontinuities in the interplanetary 
magnetic field.- One such event, recorded when Explorer 35 
was near the axis of the cylinder, is shown in Figure 2.4*. 


* 

The data in Figure 2.4 and the calculation in Figure AP-III 
are presented in a right-handed geocentric solar ecliptic 
coordinate system with +X toward the sun and +Z toward 
the ecliptic north. 
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AX so shown is the observation* of the same event by Ex- 
plorer 33 about fifteen minutes later, while this 
satellite, orbiting around the earth, was outside the 
bow shock . 

The difference in the jump in the Z components 
probably arises due to the finite 8 of the solar wind. 

An interpretation of the other two characteristic 
differences, namely the dilatation in the rise time of 
the Z component and the small peak in the X component, in terms 
of a signature of a conductive lunar interior was , : - ; - 

ruled out by Sonett et al . (197 Id), Instead they sug- 

gested that since the two satellites are widely separated 
(- 53 R„) , the signal difference might be attributed to 

iii 

a natural difference in the*- solar wind field at their * 
locations. Though such solar wind field differences might exist, 
we examine the possibility in Appendix III that these 
characteristics are caused by the surface wave nature of 
the field in the void region. We find good qualitative 
agreement with the data using values predicted by Equa- 
tion (2.4.1). 

We now return to our derivation of the field and 
proceed to complete our representation by discussing the 
"end effect" field. 
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At the vacuum plasma interface, the field must have 
a. vanishing normal component of its magnetic field and 
tangential component of its electric field* This comes 
about because we assume the solar wind to be unperturbed 
outside an infinitesimally thin boundary layer at the cylinder 
boundary. The void region is thus similar to a hollow 
pipe with perfectly conducting wall. A possible solution 
of the Maxwell equations in that geometry can be given 
in terms of the cylindrical waveguide modes (see, for 
example, Stratton, 1941). That is: 

a - 

E= 2 V x a + VxVxO a 
,o z £m z 


where 


a .a . /Sin-, , 

A *m lu,p o J m (B £m p/ V e { cos )rn * 


“ B L J m (a £m p / R M> e ’ ato 


?/RM { sin }m<|, 
cos r 


(2.4.2) 


where and a„ , 8, = 1, 2, are the roots of 

Jim 3cm 

jj^(z) and J m (z) respectively and where aJ m and 

are the coefficients of the cylindrical TE and TM modes 

respectively. 

Equation (2.4.2) is valid when the frequency is 
well below the cut-off frequency for mode propagation. 
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The smallest cut-off frequency for either the TE or 
C ^ 11 

TM modes is ~ ^ 5 0 Hz, where c is the speed of 
2itR m 

light. Since we are interested in frequencies much lower 
than this, the approximation is justified. 

The mode representation given by Equation (2.4.2) should 
be sufficiently accurate downstream of the anti solar point 
but extension of its validity to the lunar surface involves the 
so-called Pay leigh hypothesis (Millar, 1973). Indeed, in the 
vacuum region bounded by the lunar ■ surface and a plane perpendi- 
cular to the axis of the cylinder and containing the antisolar 
point there is no physical reason to impose the condition that the 
modes are decaying downstream. To illustrate this point, let us assum 
there is a magnetic dipole on the axis of the cylinder half-way 
between the center of the moon and the antisolar point. 
Furthermore, let us assume that the sphere is inside an 
infinite cylinder. Solution of this problem involves 
the Green function associated with a magnetic dipole 
inside a cylinder (see, for example, Smythe, 1968). This 
Green function is expressed in cylindrical TE and TM 
modes that are decaying upstream for z < R^/2 and down- 
stre, a for z > R^/2 . Though, in that case, the field 
can be represented by both upstream and downstream 
decaying modes it does not necessarily follow that in 
the region z > 0 and (z + p ) / > we cannot ex- 

press the field with only nodes decaying downstream. Indeed, 
this is usually not the case. We recall that the poten- 
tial field outside a sphere of radius R,,, caused bv a 
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distribution of charge and currents inside it can always 
be represented by a series of electric and magnetic 
multipoles situated at its center. For example, for 
the problem above, we can find an equivalent distribution 
of magnetic multipoles at the center of the sphere, and 
then use the Green function associated ■ with each multi- 
pale to represent the total field outside the sphere. 

Since the equivalent multipoles are situated at the 
center of the sphere, only +z decaying modes are re- 
quired in 2 > 0; but, on the other hand, the expression 

2 2 1/2 

is valid in general only in (z + p ) > 

Let us formalize this concept for the case of the 
moon. We assume that the field in the downstream cylinder 
caused by the currents and charges inside the moon and 
on its sunlit surface can be represented by a series of 
electric and magnetic multipoles situated at its center. 
This partial field is the elementary solution of Max- 
well's equations and can be found in most text books 
on. electromagnetic theory. To the field of each multipole 
we can add a series of spherical TE and TM modes, regular 
at the origin and such that the total field at the 
boundary of the cylinder has a vanishing tangential 
electric field component and normal magnetic field 
component. This field can be found through the use of 
a dual Fourier -Bessel series, as shown in Appendix IV. 

The result, in a spherical coordinate system, is as 
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follows 


H = 


V x E 
iwu 


E = 


.2 . a 


Z lajy R., t [V a + V x V x 0 a] 

_ H o m nm y nm r nm r 

n ,111 
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where the constant N m (2&) and M m (2£) are defined in 
Appendix V. 

Following Watson (1930) , we can show that the 
potentials ^^and are absolutely convergent pro- 

vided r < 25^. If we are interested in the field farther 
downstream, we can transform our. multipole representa- 
tion into a series of cylindrical waveguide modes. The 
algebra necessary to accomplish this transformation is 
outlined in Appendix IV and the result is as follows : 


± V x E 


H = 


E = 




V x + V x V x + V x V x A^a. 

nm z nin z nm z 


(2.4.4) 


where 
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B. 


ml 


= I 


2t° 

nm m m.2/ 


n (n-m)! 


n j (J / v n~2 

. _ . 2nR M d nm < “na ) 

Vi i / N , T 2 , X 

n (n-m)! V^a^) 

where t a and d a are the multipole coefficients given 

nm nm 

by Equation (2.4.3) and which must be found by applying 
the boundary conditions at the lunar surface. 

We should . point out that in this expression, 
a formal interchange of the order of summation between n 
and l was made. However, for a given 

n, if the l summation is made first, the cylindrical wave 
guide mode expression for the field of a given multipole 
is poorly convergent for small values of z and, in fact, 
usually diverges for z = 0. Thus it is difficult to 
justify theoretically the 

use of Equation (2.4.2) to express the end effect field 
in the void. This difficulty can, however, be countered 
by checking how well our boundary conditions are matched 
on the back side of the moon and, if' a good match occurs, 
appeal to the uniqueness of the solution for Maxwell's 
equations to ascertain the validity of the representation. 
Both the^representation given by Equation (2.4.2) and 
Equation (2.4.3) were used to match the boundary condi- 
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tions behind the moon and in both cases a good fit was 
found. 

Schwartz et al. (1973) also used the cylindrical 
waveguide mode to represent the "end effect” field and 
they also found a good match in the boundary conditions 
(private communication from Schubert). 

2.5 The Boundary Conditions on _ the Lunar Surface^ 

On the downstream side of the moon, all the com- 
ponents of the magnetic field and the two tangential 
components of the electric field must be continuous at the 
void/moon boundary. In addition, the normal component 
of the electric field must be zero just inside the moon 
since we assume the conductivity of the void to be 
exactly zero. But it should be noted that inside the 
void, on the lunar surface, a normal component of the 
electric field can exist due to a distribution of electric 
charge on the surface of the moon. 

On the upstream hemisphere of the moon, the tangential 
electric and normal magnetic field components must be, as 
usual, continuous. Moreover, in the frequency range 
where lunar induction has been measured the hydromagnetic 
disturbance cannot propagate upstream in the solar wind. Thus, 
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the above components of the field must assume essentially 
their unperturbed incident solar wind values on the sunlit 
surface of the moon. A partial check of these latter 
boundary conditions can .be inferred from Figures 1.2 and 1.4 
where we note that the amplification of the normal mag- 
netic field component is indeed nearly equal to one 
on the upstream hemisphere of the moon. Another partial 
test .of these boundary conditions is provided by comparing 
the data obtained from the Apollo 15 subsatellite to the 
one of Explorer 35 when both instruments were on the up- 
stream side of the moon (Schubert , 19 74a) . Inherent in the 
above boundary conditions is the assumption that a surface 
current exists within a thin boundary layer above the lunar 
surface-. which shield the solar wind from any up- 
stream influence of the moon. The Apollo 15 subsatellite 
magnetometer measurement , at an altitude of 100 km above the lunar 
surface did track rather well the unperturbed 
magnetic field observed at Explorer 35 , but the tracking 
was not perfect. The level of 

the high frequency fluctuations (f £ 0.01 Hz) was often 
seen to be somewhat higher at the subsatellite than at 
Explorer 35. However, differences in the characteristics 
of each magnetometer could also account for some of the 
differences in the magnetic records.' We should 
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noted that an imperfect confinement of the induced field 
can also be surmised from the tendency, of the amplification 
of the normal component of the magnetic field to fall 
slightly below unity in the high frequency range (Figure 1.2). 

The tracking of the field measured by the two satellites was 
often very poor when the subsatellite was near the ter- 
minator. But/ confinement near the terminator cannot be 
expected to be perfect since the dynamic pressure of 
the solar wind is nearly tangential to the surface in 
that region. We do not know yet to what extent the concept 
of perfect confinement should be relaxed in order 

to account for these observations. Nevertheless, 
this concept provides a good first-order ap- 
proximation of the real boundary conditions. 
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, Chapter III 

N umerical Solution 
Methods and Results 

3.1 Introduction 

The theory exposed in the last chapter can be ap- 
plied whether or not the toroidal H field is an important 
contributor to the total magnetic field. However, when 
the toroidal H magnetic field can be neglected compared 
to the poloidal H magnetic field, the evaluation of the 
magnetic induction is simplified substantially. In 
the next section, we examine the feasibility of ignoring 
the toroidal field component in our calculations. We 
then proceed to describe in detail the numerical method 
we used to compute the magnetic field and to discuss 

some of the numerical error resulting from the truncation 
of our various series representations of the field. In 
the last section, some of the main characteristics of the 
computed response, as a function of the parameters of 
the source field, conductivity models and LSM position, 
are examined in some detail and compared with the data. 

3 . 2 The Boundary Conditions and the Toroidal H Field 

In Section 2.5, we mentioned six boundary con- 

ditions that the electromagnetic field must satisfy on 
the dark side of the moon and three more on the upstream 
side. However, classical theorems in electromagnetic 
theory (e.g., Mtiller, 1969) show that it is sufficient 
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to impose the continuity of the two tangential com- 
ponents of the magnetic and electric field on the down- 
stream side of the moon together with the continuity 
of the tangential E components on the sunlit side in order to fully 
determine the electric and magnetic field in each re- 
gion of the lunar environment. Instead of embarking into 
a program that tries to systematically use these boundary 
conditions to compute the field, we examine their pro- 
perties in some detail since they turned out not to be 
very convenient for the problem at hand. This will lead 
us to a more proper set of boundary conditions together 
with a somewhat simpler method of solution. 

We first note that, due to the asymmetric plasma 
environment of the moon, we cannot extract from the 
boundary conditions at the lunar surface a subset pertaining to 
each electromagnetic mode, as was the case at the boundary of each 
internal layer inside the moon. Instead, they provide 
a coupling between these modes. However, this coupling 
might turn out to be poorly represented if the continuity 
of the tangential E and H components are used to find 
a numerical solution of the field. To illustrate this 
point, let us consider the continuity of the two tangen- 
tial components of the magnetic field on the downstream 
side of the moon. In order to simplify the algebra, let 


us assume that the incident field has its normal parallel 

SSSiSSs 
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to the axis of the' cylinder ( ip = 0 and say H qx = O', 

H = 1 in Equation (2.2.1) . For 0 ^ 9 ^ tt/2 and supressing 
o 

the corrmon $ dependence, these boundary conditions can be written 
as follows : 


a) Continuity of H 
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In these equations we have used the cylindrical waveguide modes • 
representation of the end effect field but the multipole 
expansion would not alter the following argument. 

We note that if these two conditions were met exactly, 
we could operate on both sides of the first equation with 

. and subtract the second equation to obtain 

9 0 

For 0 -< 9 tt-,/2 

(n) (n+l)T pJ(cos9) 

X sin 6 = 0 

n ik // 

This last equaion is equivalent to imposing 
the constraint that the normal component of the electric field 
associate with the toroidal H mode is zero on the dark 
side of the moon. Thus we infer that the continuity of the 
two tangential magnetic field components inherently 
implies this condition. However, in practice, we can 
only satisfy these boundary conditions approximately, 
so the validity of this inference needs to be investiga- 
ted from the numerical point of view. Let us consider 
the order of magnitude of the various terms in our two 
equations. The poloidal H and cylindrical TE field are 
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either of the order of the incident field or of the 
order of the toroidal H field, depending on the one that 
dominates. An estimate of the toroidal field can be based on 
the data of Figure 1.2 which shows that the solar wind 
magnetic field is perturbed only slightly on the lunar 
surface at low frequencies. Because the amplitude of the 
toroidal H response should be maximum at low frequencies 
(see, for example, sill , 1970), we infer that this field 
must be much smaller than the incident field over the 
whole frequency range. This last inference is given strength 
by the fact that the low frequency disturbance can also 
be attributed to either a small poloidal H signal or to 
noise which seems to indicate very 
low near-surf ace conductivity* . if the 


* 

In our discussion, we have assumed that the relaxation 

time associated with charge diffusion e /a) is much 

smaller than the period of interest. For example, at 

-14 

T=25 sec, we must have a>>4x!0 “ mho/m. Strangway (1968) 

has suggested, however, that near surface conductivity on 

— 1 3 — 1 6 

the moon might be as low as 10 -10 mho/m. in that 

case displacement current would dominate over the conduction current and 
not only the expression for the toroidal H field would have 
td'^b'e revised but also its associated boundary conditions. 
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toroidal H magnetic field is very small compared to the 
poloidal magnetic field we can not devise easily a numer- 
ical algorithm to match our four boundary conditions and at 
the same time conserve their independence. However, if 
we are interested in calculating only the magnetic field 
and are willing to accept an error in doing so of the order 
of the toroidal H magnetic field, we can simply drop the 
toroidal contribution from the continuity condition on the 
tangential magnetic field components. Then, using the 
tangential components in conjunction with the boundary 
conditions on the normal component of the magnetic field 
we can determine the magnetic field completely. This 
then seems to be the proper set of boundary conditions to 
use since solving simultaneously for the electric and 
magnetic field would only provide a correction of the 
order of the neglected toroidal H magnetic field which 
is likely to be small if not negligible. 

3 - 3 Numerical Method and Precision of the Solution 

As mentioned in the preceding section, the magnetic 
field can be completely described by matching only the 
boundary conditions associated with its tangential and 
normal components, neglecting terms involving the toiroidal 
H field. in order to obtain this solution numerically, 
we note first that, for each value of 0, each component 
of the field can be 


-71- 


expressed in a Fourier series in the <p coordinate. The 
Equations (2.3.1), (2.4.1), and (2.4.2 or 2.4.3) are 

already in that form and we show in Appendix VI how the 
normal component of the incident solar wind ^magnetic field 
can also be expanded in such a series. Because of the 
orthogonality of these different harmonics, the partial 
field associated with a given pair (m,a) in one region 
of the lunar environment needed only to be matched with 
the partial field associated with the same pair in another 
region. Consequently, the <J> dependence can be suppressed 
from our boundary condition and we are left with the 
task of matching only their polar angle dependence, which 
we can do by writing down our boundary conditions at K 
different values of 0 both on the sunlit and dark side of 
the moon and then use . only a finite number of parameters 
to represent the poloidal H and "end effect" field as- 
sociated with the pair (m,a) . The resulting system of 
equations can then be solved by the method of least 
squares. This process can be summarized by the following 
matrix equation: 


C X = S 
mo raa me 

where , for reasons to be discussed shortly, 
are distributed as follows: 

* 


(3.3.1) 
the elements 
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1. The top J elements of the column of parameters 

X are filled with the unknown coefficients of the end 
me 

effect field (or t^ m ) and are followed by the co- 

efficients of the poloidal H field (a^) . 

2. The top 3K lines of are used to represent 

\u\j 

the difference between the end effect and poloidal H field 
of the three magnetic field components on the dark side 
of the lunar surface whereas we reserve the last K lines 
to represent the normal component of the poloidal H mag- 
netic field on the sunlit side. We note that, though in 
theory we need to use only one component of tangential 
magnetic field, we found that in practice, a more uni- 
form distribution of the truncation errors result when 
both components are used. (The case m = 0 is somewhat 
special since there is no component associated with 

this value of m, so only the continuity of H and H are 

o r 

used in this case) . 

3. Our boundary conditions require that the top 3K 

elements of the column S be filled by the matching 

me 

values of the surface wave field whereas its last K elements 
are reserved to express the normal component of the in- 
cident solar wind magnetic field at the same K values of 
0 used on the left hand side. 
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By using a series*- of Householder . transformation to 
solve Equation (3.3.1) by the method of least squares 
(see, for example, Golub, 1965), we realize a substantial 
economy in the computation when, we exploit the following 
three properties of the elements of the first J column 
of C ma (associated with the end effect field) . 

a) Each of these elements are real as compared to 
the generally complex value of the elements associated .with 
the poloidal H mode. 

b) The last K elements of each of these columns are 
filled with zeros since the sunlit side magnetic field 
does not involve the end effect field. 

c) These elements are independent of both the fre- 
quency and the conductivity model used [as can be seen from 
Equation (2.4.2) or (2.4.3)]. 

This last property is especially useful when we need 
to evaluate the field at different frequencies and for dif- 
ferent conductivity models since the Householder trans- 
formation associated with the first K columns may be saved 
and then used repeatedly in each of these cases. 

Before investigating some properties of the field, 
we need to specify the numerical error associated with the 
various truncations in our series representation of the 
field and by our choice of a finite number of points to 
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match its boundary conditions. With regard to this last 
point, we found that the field obtained by matching the 
boundary conditions at five degree intervals in the 6 co- 
ordinate did not differ by more than 1% from the one 
computed if three degree intervals were used. Thus, the 
former choice was judged sufficient for our purposes. We 
therefore restrict our discussion to the error arising from 
the following two main types of truncation: 

1. A truncation in the number of .Fourier harmonics 
(m = 0 , 1 , . . . , M) . We discuss the criteria .used to 
effect this truncation in Appendix VI and show that a 
relative error of less than a few percent should result. 

2. For each Fourier harmonic specified by a pair 
(m,a) we need to truncate the series representations of the 
poloidal H and end effect field. 

Before proceeding to discuss this last point, we 
should mention that various other types of truncation 
are generally required to evaluate each element of the 

matrices C . and S _ (e.g., truncation in the series 

representation of a given function). Their associated 
error, however, can be rendered negligible compared to 
the ones above by using sufficiently accurate algorithms. 

Let us define the following two measures of the 
relative mismatch in the boundary conditions arising from 
our choice of a finite number of unknowns to represent the 
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poloidal H and "end effect" field. 


R i?(e) = 


R 2 i < 9 1 = 


mo , mo 
a i - b i 


\ma 

a i 



ma 


, ma 


a i 


b . 

i 


ma 
a . 

1 


where a^ a (0) = the i component of the poloidal H field 

on the lunar surface (i — r, 0, 4>) as- 
sociated with the (m, o) harmonic 
b™ a (0) =1) On the dark side: the i component of 

the surface wave plus end effect field on 
the lunar surface associated with the (m,a) 
harmonic 

2) On the sunlit side: the i component of 
the incident solar wind magnetic field 
on the lunar surface associated with the 
(m,c) harmonic. 

Using the triangle inequality, we can show easily 
that 

To illustrate the re suiting mismatch in the boundary, condition 
we shall concentrate on the case when the incident field has its normal 
parallel to’ the .axis of the cylinder (and, say, H = 0) . This 

U A 


R^(0) * RJl ( 6 ) . 
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source field involves only the pair (m,c) = (1,1), but 
the results are typical of those obtained for other 
values of m. We shall show shortly that, for a given 
number of field parameters to represent the field, the 
mismatch tends to increase with frequency. Thus,, we shall 
examine in detail only the high frequency case, f = 0.04 Hz, thereby 
setting an upper bound to the error we might expect. The mis- 
match also depends to some extent on the conductivity 
model used. The results discussed below are typical of 
models that give good fits to the data but other 
models used in the text generally will agree with this 
error estimate to within a factor of 2. 

We illustrate in Figure 3.1 the result obtained when 
the poloidal H and end effect field are represented by 
10 unknowns each. The component on the dark side is 
shown together with the two measures of the mismatch 

defined above (dot for log R^) . We note that either 

2 <P 

measure is indicative of a relative error of less than 1%. The 
mismatch in the component is examined in Figure 3.2. 
This time, however, both measures of the mismatch are 
log but the dot displays the relative error obtained 

when 15 parameters are used to represent each field (in- 
stead of 10 for the square) . Though this increase in the 
number of parameters reduces the relative error by about 
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a factor of 2, the 1% accuracy obtained with the smaller 
number of parameters was considered sufficient for our 
purpose. 

We have not distinguished above which of our two "end 
effect" field representation was used. The reason for* 
this is made plain in Figure 3.3 where we illustrate the 
results obtained from both expressions. In each case, ten 
parameters were used to express the "end effect" and 

poloidal H field and the measure of the mismatch shown 

11 

1S R 2r * T ^ e re l atlve error using the multipole ex- 
pression [Equation (2.4.3)] is represented by the dots 
whereas the triangles exhibit the one obtained when the 
cylindrical waveguides mode [Equation (2.4.2)] are used. 

Not only is the error in both cases nearly identical', 
but also the computed field does not differ by more than 
1% (approximately the radius of the dot in the upper 
figure). Moreover, similar results were obtained for the 
other components of the field. These results suggest two 
conclusions. First, though in our argument in the last 
chapter we suggest that the multipole repre- 
sentation is theoretically more adequate for the problem at 
hand, the numerical results do not support this contention 
Indeed, the closeness of the computed field using either 
mode of expression hinted that they are merely equivalent 
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Figure 3.3 
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ways of expressing the field and that the Rayleigh as- 
sumption embodied in the cylindrical waveguide mode 
expression is justified. Secondly, the closeness of the 
relative error to each representation in addition to the 
fact that their amplitudes are approximately the same on j- 
the sunlit and dark side of the moon suggest that the 
truncation j.n the poloidal H expansion, is the more important 
contributor to the mismatch. In other words, for a given 
number of parameters , the error can probably be reduced 
by using more parameters in the poloidal H field repre- 
sentation than in the "end effect" field representation. 

This could have been anticipated on the basis of the im- 
portant role the day/night asymmetry ought to play in the 
computed response. 

Since the 5% relative error in the radial component 
relative error in the radial component shown in Figure 3 . 3 
has little consequence on our future results, all the 
calculation in this work • (except one to be mentioned in 
the next section) were done using ten parameters to 
represent the end effect and poloidal H field. We should 

«v -4 

keep in mind that the relative error illustrated in the 
preceding three figures are upper bounds on what we can 
expect in the frequency range of interest. We show in 
Table" 3.1 typical results obtained at lower frequencies for 
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the same number of parameters as quoted above 


Table 3.1 
- Relative Error 

11 „H s R 11 

f (hertz) K 2 <|> % R 26 % _2r 


0.0002 

0.005 

0.05 

0.05 

0.002 

0.05 

0.5 

0.7 

0.02 

0.1 

0.8 

3 

Though 

we can deduce 

from this table that, 

at low frequencies 


an acceptable precision can be reached with a substantially 
lesser number of parameters to represent the field, we did 
not use this option in the work. 

3 . 4 Propert i es of the Field 

The data presented in the first chapter (and in 
Section 2.3), depend not only on the internal electrical 

conductivity of the moon but also on: 

1. The detailed characteristics of the source field, 

2. The exact position of the magnetometer measuring 

the field on the lunar surface, 

3 . The presence of noise associated with processes 

not accounted for by the theory. 

Unfortunately, the data are averages over an 
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unknown distribution of these parameters. In addition, 

• • VI- - 

an empirical correction was applied on the sunlit side "of 
Apollo 12 data to account for what appears to be a noise 
associated with the presence of a substantial remanent 
at this site although a theoretical foundation for this 
correction is still lacking. Therefore we proceed to 
illustrate to what extent the measured and computed res- 
ponse are sensitive to these various parameters and noises 
in order to gain insight into the correctness and limitation 
of the assumptions we shall use in the inversion process 

Figure 3.4: We show here the A . value discussed 

in Section 2.3 and A data (from Figure 1.2) foth of which 
were collected at the Apollo 12 site when the LSM was on 
the sunlit side of the moon. Also shown is the computed 
field at the subsolar point using the symmetric plasma 
assumption (dot) and our more accurate calculation in- 
volving the asymmetric plasma environment (^). 

In both cases the three-layer model given below was 
used. 


R adius (km) 

0 - 1206 
1206 - 1527 
1527 - surface 


c (mho/m) 
1.9 x 10“ 
5.12 x 10 
1.0 x 10“ 


This model was obtained by Sonett et al . (1972) as 
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a member of a set giving a good fit to the A . value when 

mm 

using the symmetric plasma assumption to invert the data. 
The source field used is impinging normal at the sub- 
solar point (i.e. f $ - 0) , and the solar wind velocity 
used was 400 km/sec. 

We note that the inclusion of the void region behind 
the moon substantially lessens the predicted lunar response 
on the sunlit side, especially at the shorter period. This 
effect occurs mainly because the confinement current on the 
boundary of the cylinder is less efficient than the one in- 
herent in the symmetric plasma model in amplifying the 
front side response. In order to build up the response to 
the values, a higher conductivity distribution is re- 

quired near the lunar surface. However, we also note that the 
extraction of the values from the sunlit side data involve a 

relative correction much larger than the correction 
associated with the introduction of the void region behind 
the moon. Thus, though our theoretical refinement does 
permit unification in a single framework of both the sunlit 
and dark side responses, its ultimate superiority over 
the symmetric plasma theory in obtaining an accurate lunar 
conductivity model is largely dependent on the accuracy 

of the correction made to obtain the A . values 

mm 




DEPTH ( km) 

Figure 3.5: This is a plot of three models 

proposed by Ward (1969) as a guide for the range of con- 
ductivity .-'expected inside the moon. He bases his 

estimate of the model labelled "B" on a lunar temperature 
model proposed by Fricker et al . (1967) and assigns a 
conductivity at a given depth consistent with, the value realized in 
the earth's mantle at the temperature characteristic of 
this depth. Due to various uncertainties as- 
sociated with this correlation. Ward proposed models A and 
C as reasonable error bounds. Actually 

these bounds are wide enough to embrace not only all the 
models giving a good fit to the electromagnetic data but 
also practically all the models based on the various lunar 
temperature distribution, composition, and conductivity- 
temperature relationship proposed in the literature. Thus/ 
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Ward's models formed a good starting point to investigate the 
sensitivity of the data to the conductivity distribution. 


& 



Figure 3.6: We show the subsolar point response of 

the preceding three models as a function of frequency, when 
the source field parameters are ^ = 0 and = ^300 km/sec. 

Also plotted are the A . values*. 

■ L min 


* 

In all the following plots of amplif ication vs. frequency, 
the frequency scale is logarithmic and the response was 
computed at the frequencies specified by the ticks on the 
bottom scale. These are: 0.0002, 0.0005, 0.001, 0.002, 

0.005, 0.01, 0.02, 0.04 Hz . In some graphs,. the lowest 
frequency is omitted. 
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The low frequency response is indicative of a deep 
lunar conductivity similar to that of model C, 
whereas the high frequency response is indicative of a conduc- 
tivity structure at shallow depths similar to that of 

model B. Thus, we infer a moon of fairly 

. . —4 -2 

homogeneous conductivity between about 10 to 10 mho/m 

below approximately 100 km covered by a fairly resistive 

layer. A more explicit description of the conductivity 

structure will be found in the next chapter. 



F RE CX UENC Y (hertz) 


Figure 3.7: The same source parameters and three 

conductivity models discussed above are used here to compute 
the antisolar point tangential response. 
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Clearly, the large errors associated with these data 
and their limited low frequency coverage range 
render them far less useful than the more extensive 
sunlit side data in discriminating between conduc- 

tivity models. To make matters even worse, as we shall see 
shortly that the high frequency data (above 
approximately 0.005 Hz) is very sensitive to the parameters 
of the source field and LSM position. Nevertheless, the 
two data points at the lowest frequency might be symptomatic 
of a difficulty associated either with the theory or a 
source of noise unaccounted for by it. We note that, from 

this and the preceding figure, the computed antisolar point 

/ 

response is distinctly smaller than the one at the sub- 
solar point (by a factor of about 1.6 at this frequency 
for model B) . Yet the two data points at 0.0017 Hz do 
not seem to bear this out. Indeed, the A^ value is only 

about 1.1 times smaller than the A . value at this fre- 

mm 

quenc'y, whereas the A value is actually 1.1 times greater. 
However, due to the large uncertainties associated with 
these measurements, we w refrain from drawing definitive con- 
elusions as to their consequences. 
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ductivity structure. The source field was again speci- 
fied by = 0 and = 300 km/sec and the three con- 

ductivity models used are detailed in the figure. We 
note that the model exhibiting a conductivity 
region to about nine-tenths of the lunar radius (158 0 km) 

agrees fairly well with the A^ ±n values. Actually, this 
model was proposed by "KuckeS (1971) in an attempt to fit 

the sunlit side A and A v data using the symmetric plasma 

z y 

theory. Although the asymmetric theory lessens the 

■^•-computed response obtained through the symmetric theory, 

the A . data are also lesser than the A and A data, 
min c i 
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thus leading in each case to a similar conductivity model. 



FREQ. UENCY (hertz) 


Figure 3.9: ‘ This is a plot of the antisolar point .. 
response for the source parameters and model described in 
the preceding figure. The large uncertainties associated 
with the measurements render them practically useless 
to discriminate between these three different conductivity 
models. 

Before going any further, let me interject a few 
general remarks. 

Since the LSM at the Apollo 12 site is situated at 
2.95° S latitude in selenographic coordinates, 
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negligible errors result if we assume it to be in the 

equatorial plane. An immediate consequence of this location 

is that the north-south, east-west, and vertical data 

values (A ,' A , A ) correspond respectively to the cora- 
z y x 

puted A,, A_ , and A_ amplification factors. However, it 

<p 0 R 

should be pointed out that since the computed A^ , A 0 , and 

A*, can depend on the angle the equality will be met, in 
R 

general, only when the wave normal is in the equatorial 
plane. We did not differentiate in the preceding plot 
which angle (p and which component A Q or was used. 

This is because when the incident field has its normal 
parallel' to the axis of the cylinder, one can easily show 
that the amplification is isotropic at the subsolar and 
antisolar point. As a matter of fact, we can also show 
that, for this source field, the three amplification fac- 
tors are independent of the angle for a general value of 
0, though the two tangential amplification factors need not be 
isotropic in that case. 

In the remaining graphs in this chapter, we shall try 
to illustrate some aspects of the dependence of the com- 
puted response on the parameters of the source and LSM 
position. Though these results also depend on the par- 
ticular conductivity models used, a very good general idea 
can be gained by considering only one structure. In 





order that the results remain relevant for 


future 


work, we have chosen the model detailed in Table 3.2. Thi 


model turns out to give our best fit to the A . value 

.* . mm 

at the subsolar point, when the source field is specified 
by 0 =0 and V gw = 300 km/sec. 

Table 3 .2 

Conductivity Model 
Lunar Radius (km) a (mho/m) 


1020 

0.11 

X 

10- 2 

122 0 

0.81 

X 

io “ 3 

1320 

0.15 

X 

10“ 3 

14 20 

0.92 

X 

10 -3 

1520 

0.51 

X 

10 -3 

1570 . 

0.23 

X 

10~ 3 

1620 

0.11 

X 

10' 3 

1670 

0.42 

X 

10~ 4 

1730 * 

0.1 

X 

10“ 10 


1730 km was used instead of 1738 km but this 
has little consequence. 


We shall also concentrate our attention only on 
the response obtained in the plane formed by the wave 
normal and the axis of this cylinder. This distinction 
is irrelevant when 0-0 but in other cases I shall call 


this plane the plane of incidence. 


By symmetry. 


the H 


<P 


field in that plane depends only on the polarization 
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specified by H qx in Equation (2.2.1) whereas the H Q 
component depends only on the H q polarization. For con- 
vsniencG ; we shall assume both polarities to have unit 
amplitude (i.e., H q = =1). 



Figure 3.10: This is a plot of the absolute value 

of the <j> component of the magnetic field vs. the polar 
angle 9 (9=0; antisolar point), for different fre- 
quencies. The source field is specified by i' - 0 and 
V_._ - 300 km/sec. Because of our previous assumption, 
'?fehe power in the § component - of the source field is unity 
in the plane of incidence so the component is also 

directly the amplification factor. 
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The day-night asymmetry in the calculated response 'is clearly 
evident in this figure. However, we note that the trans- 
ition between the two responses is fairly sharp but occur- 
ring mainly on the dark side of the moon. Thus we are . left with a 
fairly extensive region near the subsolar point where the 
response is practically independent of position. 



FREQUENCY (hertz) 


Figure 3.11: This is a plot of the same A, as 

the preceding figure but this time as a function of 
frequency and for selected values of 8. Note that within 
45° from the subsolar point the amplification does not 
differ by more than a few percent. However, at 0.04 hertz and at 


45° from 
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the antisolar point the response can reach a value of 
25% higher than the response at the antisolar point, 
highest frequency. 



point for \p - 0 and various solar wind velocities 
A decrease in solar wind velocity tends to depress 

the high frequency response. 



FREQUENCY( hertz) 


Figure 3.14: The antisolar point response is 

illustrated here for different solar wind velocities but for an 
nngle o£ incidence specified by tjj = 0. If we compare the 
high frequency response to the response presented in Figures 3.7 
and 3.8, we note that the solar wind velocity can he a 
more important factor than the conductivity in controlling 
the high frequency 7 response. 

Receding page blank not filmed 
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Figure 3.15: This is a plot of the antisolar point 


k <t> 


response for a solar wind velocity of 300 km/sec 


and various angles of incidence. Again we note that an 
increase in ' the angle of incidence tends to produce a 
similar effect to that of a decrease in the solar wind 


velocity . 





Figure 3.16: This is a plot of the absolute 

value of the 0 component of the magnetic field vs. 
the polar angle. The incident field is characterized 
by ip - 0 and = 30 0 km/ sec . 

The peak that developed near the terminator 

(9 =90°) at high frequencies is due in part to the 
distortion of the incident field by the surface wave. 
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• 

component along the axis of the cylinder is associated 
with the surface wave even though the incident field has 
no component in that direction. This component becomes 
an important contributor to the total field at high 
frequencies especially near the boundary of the cylinder. 

The H q component of the field is not directly the 

amplification since the power of the incident field 
2 

has a sin 0 dependence on the polar angle. Consequently 
one must divide H 0 by sine to obtain the 0 amplification 
factor. But this involves a division by 0 and a resulting 
infinite amplification at the terminator. However, this 
problem is somewhat artificial since it arises solely 
because. we assume the incident field to have no power along 
a given direction. In an actual situation the incident 
field has power in every direction though its distribution 
can be highly anisotropic (see Belcher et al. for a dis- 
cussion of this point) . in order to counter this difficulty 
to some extent, we shall henceforth consider only tangential 
amplification at least 45° from the assumed direction of 
zero power in the incident field. We think this procedure 
provides us with a fairly accurate picture of the actual amplification 
though we have not proven rigorously this assertion 

Before going to the next figure, we recall, 
that when the source has its normal parallel to the axis 

DBIGINAIi PAGE IS 
m POOR QUALITY 
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of the cylinder, the amplification is isotropic at the 
subsolar and antisolar point. Consequently, no change 
results if, in Figures 3.12 and 3.14 the amplification 
factor is replaced by . 


FREQUENCY {he rt?) 


Figure 3.17: In this plot, the source field is 

specified by ^ — 45° and Vt — 300 km/sec. 


Let us consider first the curves representing the 
sunlit side response. If we compare the above A 0 response 
at the subsolar point with the subsolar point response for p - 0 
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(from Figure 3.12), we notice about a 5% increase 

in the computed response at high frequencies 
Thus, in contrast to the response, an increase in the angle 
of incidence tends to increase somewhat the subsolar 
point A a response. When = 45°, the response is aniso- 
tropic at the subsolar point with the Ag response some- 
what larger than the A^ response at high* frequencies. 

There is still a small anisotropy at 45° from the sub- 
solar point (0 = 135°) but it does not amount to more than 
about 5% though it is still A Q that generally dominates 
the response. Based on these preliminary results, 
we would expect anisotropy to be observed on the 

sunlit side with A Q generally larger than . How- 

ever, if we examine the data of Figure 1.2, we see 

that A (i.e., A^) is generally higher than A (i.e., 
z <p y 

A ) except at high frequencies. 

0 

However, at low frequencies, the remanent field 
seems to play an important role in determining the aniso- 
tropy . Thus we cannot conclude at this point that our pre- 
liminary calculations are contradicted by the observations. 

By comparing Figures ; 3, 15 and 3. 17, we note that at the 
antisolar point, A Q , like A., tends to decrease with &n 

' U(p ’•> 

increase in the angle of incidence. When \p - 45°, the 
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compute d response at the anitsolar point and at 45° from 
it has a somewhat mixed character, at the antisolar point 
A , dominates while at 45° from the antisolar point A fl is 

(p 0 

the larger of the two components. Due to the order of 

magnitude of the differences involved, however, we would 

suspect that A Q is generally the larger response for data 

taken within the first 45° from the antisolar point. On 

the other hand, the data of Figure 1.3 exhibits values of 

A larger than A , contrary to our expectations, and we 
z y 

cannot here involve the remnant field to explain this 
divergence. But since the mean square errors associated 
with the average value of the measurements overlap each 
other, these uncertainties preclude any definite conclusion 
at this stage. 



3.18: This is a plot of the radial amplification 

factor as a function of position on the dark side when 
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the source field is specified by = 0 and Vg ^ = 300 km/sec. 
The values shown start at 10° from the anti solar point 
in order to avoid the difficulty associated with the fact 
that the incident field has no power along the radial com- 
ponent at the antisolar point. At the frequency 0.04 Hz, 
we have plotted the response computed using 10 (heavy dot) 
and 15 parameters in our 'representation of the poloidal 
H and of the end effect field. Though the latter case 
smoothes out the response, the relative amplitude of the 
oscillation stiM present in the former case does not 
reach more than 5% of the total field. Note that the 
radial response is characterized by a broad minimum with 

practically constant value within 45° from the anti- 

„ ■&> * 

solar point. Thus variation of the response with position 
is not very important in that region. 
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Figure 3.19: This is a plot of the radial ampli- 

fication using a source field specified by iii = 0 but 
with different solar wind velocities. The amplification 
values shown were calculated at 25° from the antisolar point 
and are average at' high frequency to minimize the effect of the 
oscillation shown in the preceding figure. We note that 
the response decreases from unity at low frequencies, 
reaches a minima, and then starts to increase at high 
frequencies. At still higher frequencies, the response 
reaches a maximum and then decreases monotomically 
to 0 . This behavior is typical only of low angle of 
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incidence source fields and is related to the distortion 
of the incident field by the surface wave. However, 
such behavior is probably not typical of what is observed 
near the antisolar point, as shown in the next figure. 



Figure 3.20: 

various angles of 

v T = 300 km'/ sec. 
SW 


This is a plot of the A R response for 
incidence but constant solar wind velocity: 
The position of the LSM was assumed to 


3 00 km'/ sec. 
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be at 25° from the antisolar point along the curve 
$ = 270°. The computed response at any polar angle 
in the range 0 = 0 - 45° is typically within 5% of the 
radial amplification factors exhibit in Figure 3.20. 

(The only noteworthy exception is the case ip = 60° at 
0.04 Hz where the response varies from 0.13 to 0.20 in 
the polar range 6 - 0 - 45°) . 

We note that the computed response is very sensitive 
on the angle of incidence and thus the data hardly can 
be inverted without some information about the distri- 
bution of this quantity in the incident field. The fact 
that the data tend to agree with curves of the computed 
responses characterized by large angles of incidence 
probably comes about because near the antisolar point the 
power in the radial component is heavily biased toward 
such angles since they are associated with relatively 
large magnetic field components along the axes of the 
cylinder . 

Before proceeding to the inverse problem let me sum- 
marize the result of this section. 

Let us consider first the dark side data. For the 
radial amplification we. just showed that a reasonable 
fit to these data can be obtained by using a model fitting 
the A min Value on the front side and choosing relatively 
large angles of incidence for the source field. However, 
due to the large sensitivity of this response to the 
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angle of incidence, a proper inversion of the,* dark side 

♦ 

radial amplification data must await additional information 
on the source field. This need for additional information 
on the source field -parameters is also true for the tan- 
gential amplification on the back side, especially above 
0.005 Hz. A proper interpretation of these data is 
further compounded by our ignorance of the exact position 
of the LSM when these data were taken. For frequencies 
smaller than 0.005 Hz , this latter dependence seems to be 
the most important parameter affecting the response (apart 
from the conductivity model) . We noted that though the 
radial data and high frequency tangential dark side data 
tend to confirm the soundness of the asymmetric theory, 

& __ the same cannot be said of the low frequency tangential 
dark side data. The measured tangential response on the 
dark side is almost equal to the one measured on the 
sunlit side whereas the theory predicts that it should be 
substantially less than the sunlit side response. We 
might also point out that additional data on the depen- 
dence of the response with position were published 
recently by Smith et al . (1973). But there again, the 

measurement at 0.0017 Hz did not exhibit clearly a smaller 
response on the dark side than on the sunlit side. However, 
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th e large uncertainties associated with these data 
preclude any definite conclusion at this stage. Need- 
less to say, complete confidence in the asymmetric theory 
can only be gained by its confirmation over the whole fre- 
quency range of interest. Therefore we feel that addi- 
tional measurements of the dark side response should be 
carried on, especially in the low frequency range and 
that noise possibly arising from drift current at the 
plasma vacuum interface should be accurately- assessed 
in the future . 

The front side data are indicative of a much smaller 

relative dependence on the source parameters and 

LSM positions and also a much greater sensitivity to the 

conductivity model, especially at shallow depths. A 

major drawback, however, is the noise associated with the 

remanent field at the Apollo 12 site. V7e correct empirically 

for this source of noise bv calculating the A . values 

mm 

but this process also removes some of the natural aniso- 
tropy we might expect between the two tangential responses. 
However, the fact that the Aoollo 12 A . values agree 
with the Apollo 15 data adds to our confidence in their 
use . 


nRlGINAD PAGE ® 
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Chapter IV 

Inverse Problem and Conclusion 

* ' .i ' i 

4 . 1 Introduction 

A substantial effort to remove some of the anthro- 
pomorphic character associated with the inversion of 
geophysical data recently has appeared in the literature. 

In the next section, we shall describe the formalism we 
used for our problem. Our analysis followed 
substantially the treatment given in Wiggins (1973) . 

Two studies of the resolution of the lunar conduc- 
tivity structure provided by the magnetic data have already 
appeared in the literature (Hobbs, 1973 and Phillips, 

1972). However, they both assume a symmetric plasma 
theory and an infinite value for the solar wind velocity. 
Moreover, since these studies appeared, the frequency 
range of the available data has been extended and the 
effect of the remanent field has at least been empirically 
removed. Thus we shall endeavor hereto bring the sub- 
ject up to date. 

The third section shall be devoted to a description 
of the results obtained when the lunar conductivity struc- 
ture is modelled as a conventional layered sphere. In 
the fourth section, we shall constrain our conductivity 
model by assuming a given temperature distribution in- 

Paoe m 
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side the moon. The object then will be to understand 
to what extent the conductivity model constrained the 
conductivity-temperature relationship of the material 
likely to form the bulk of the moon. In the last sec- 
tion, we shall summarize our results, pointing out 
various directions that can be followed to improve them, 
and examine some of the geophysical constraints the 
results seem to imply. 

4 . 2 Specialized Solution of the Forw ard and Inverse 
Problem 

As we pointed out in the last chapter we shall try 

to match the data by assuming the magnetometer to be at 

the subsolar point and the incident field to have its 

normal parallel to the axis of the cylinder. In that 

case, only the harmonics specified by m = 1 are involved 

in the computation of the field. Moreover, since the 

response is isotropic at the subsolar point for the 

assumed source field, we can conveniently specialize its 

characteristic to the case H - 0 and H - 1 [see, for 

ox o 

example, Equation (2.2.2)], which in turn involves only 
the harmonics specified by a - 1. The theoretical tan- 
gential amplification factor is then simply given by 

A - ( % H// 2 | cj) = 0°, 0 = 180° (4.2.1) 

where the bar stands for the complex conjugate and 


where 
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H 


^ I 0=0 0 , 0=180 ° 


a\L P 1 (cos 6) 

£ n l n n 
n 


sinB 


6=180° 


From Equation (3.3.1), the coefficient are 

found by solving the matrix equation 


c n x ii s n 


(4.2.2) 


It is interesting to point out that even for 
an arbitrary incident field, the tangential amplifica- 
tion factor at the subsolar and antisolar points involves 
only the evaluation of the partial field associated with 
the harmonics specified by m- = l'and c — 1,2. This comes 
about because, in our chosen coordinate system, the 
partial field* associated with each Fourier harmonic in 
the <j> coordinate can be solved separately and also be - 

cause the terms of the form ~ ^ c . OS — and 

sm 6 


3P 


m 


n (cos 9) . 


30 


m our general expression of the lunar surface mag- 


netic field [see Equation (2.3.1)] are equal to zero at 
0=0° and 180° if m ^ 1. A similar result holds for 
the radial amplification factor at the subsolar and 
antisolar point but this time, only the harmonic speci- 
fied by m = 0 needs to be considered. 
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For most of the following results, a search was 
made for a parameterization of the conductivity structure 
yielding a good fit to the response at eight frequencies 
approximately equidistant on a logarithmic scale. These 
were 0.0002, 0.0005, 0.001, 0.002, 0.005, 0.01, 0.02 


and 0.04 Hz, corresponding to periods of 25 sec to about 
1.4 hour. Also, eight layers with fixed thicknesses gen- 

orally were used to model the conductivity structure. 

A logarithmic coxxectiun to an initial parameteriza- 


tion can be found through the matrix equation 


By — P 


(4.2.4) 


where 

31nA (f . ) 

O = 1 

ij 31n o . 

D 

y. = Ain a. 

r l 

and 

A ( f . ) measured 
^i ^ A (f . ) computed^ 

From Equation (4 . 2 . 2) , we note that there are two main 
steps involved in calculating the partial derivatives 
implied in Equation {4 . 2'. 4) . First, the partial deriva- 
tive of the L factors with respect to each conductivity 
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parameter must be evaluated. An analytical expression 
for these derivatives can be worked out easily from our 
expression for the L r in Appendix IX, and the numerical 
computation can be most conveniently performed together 
with the evaluation of the L n factors. We also require 
the partial derivative of each coefficient a^ with res- 
pect to each conductivity parameter. These derivatives 
can be obtained by taking the partial derivative of 
Equation (4.2.3) and noting that S-^ is independent 
of the conductivity, we obtain: 

„ 3X n _ 3C u „ 

C 11 3a i " 3a. X 11 (4.2.5) 

Once Equation (4.2.3) is solved and the partial derivative 
of the L r factors are computed, the right-hand side of 
Equation (4.2.5) is completely determined. Moreover, 
since the left-hand side of Equations (4.2.3) and (4.2.5) 
involves the same matrix C^, the set of Householder 
transformations used to find the least square solution 
of Equation (4.2.3) can be conserved and then used 
repeatedly to solve Equation (4.2.5) for the partial 
derivative terms. 

In order to gradually reduce the effect of small 
eigen-values associated with Equation (4.2.4), the 
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generalized inverse operator * 

T 2 -1 T 

(B B + e l) ET 

was used to find a solution (Maddfen, class notes) . A 
practical way of finding the solution through this opera- 
tor is to follow the method suggested by Golub (1§65) 
which consists of applying a series of Householder trans- 
formations to the modified matrix equation 


- 


— “ 

B 

y = 

P 

el 


9 9 9 

0 

— - 




The choice of e was made by a trial and error 
search for a number giving a stable iteration and 
yielding models giving a good fit to the data. We 
found that, if this value was chosen to be between 
0.05 and 0.1, which usually corresponds to allowing 
from 2 to 3 linear combinations of parameters to par- 
ticipate in the solution, both a good fit and a stable 
iteration were realized. 

Once Equation (4.2.6) is solved, the new model 
parameter can be found through 


new 
°i = a i 


old ^i 


e 


(4.2.7) 
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The probable presence of systematic errors in the 
data ruled out a rigorous discussion of the trade-off 
between the resolution of the parameter correction and 
the standard deviations associated with uncertainties 
in the data. However, since a quantitative basis for 
a discussion of parameter resolution is most naturally 
provided by this trade-off, we shall assume the ob- 
servational error at each frequency to be uncorrelated 
and to have a relative standard deviation of 5%. This 
value* corresponds roughly to the standard deviation at- 
tributed by Sonett (1974) to the A z and A^ values result- 
ing from his data analysis. However, this estimate 
seems clearly optimistic in view of the relatively large 
correction** that was applied to these data to account for 
the presence of the remnant field (see Figure 3.4). To 
alleviate some of this uncertainty, we shall use several 

versions of the data as given by the A . , A and A 

J min z y 

values. 

It is well known (Lanczos, 1961) that the matrix 
B can be decomposed as follows: 

* 

V T (4.2.8) 


B = U 


0 
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where the elements of the diagonal matrix of eigenvalues 

can be conveniently ordered such that and 

where U and V are the matrix of eigenvectors of BB T 
T 

and B B respectively. 

The resolution matrix of the parameter correction 
associated with the k largest eigenvalues of B is 
defined as 


R = v V 
K K K 


(4.2.9) 


where the K columns of are the first k columns of 

V. 


The standard deviation of the parameter correction 
associated with the uncertainties in the data is then 
given by^ 

k 99 1/2 

(STD. DEV.) = 0.05 ( l V; /AT) (4.2.10) 

i=i 31 i 


Similarly, we can define a resolution matrix of 
the observations by 


D k * 


T 

u k u k 


(4.2.11) 


where the k columns of are the first k columns of U. 
Examination of each row 'of permits us to determine 
to what extent the information contained in a given 
observation is used to estimate the k linear combina- 
tions of the parameter correction implied by R, . 

*■ k 
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Before displaying the results, let me add a general re 
mark . 

In all the models shown in this chapter, the con- 
ductivity of the top layer was kept constant in the in- 
version routine at 10 mho/m. This has little effect 
on the resulting model when only the poloidal H mode is 
used to invert the data, but it effectively renders the 
toroidal H magnetic field negligible, in harmony with 
the data. 

4 . 3 Results from the Inversion of a Layered Sphere 

Without Thermal Constraints' 

Since we feel that the reader should get 
an overview of the results, before the significance of some 

of the details can be appreciated, we shall limit ourselves 
here to their description, leaving a discussion of their 
common characteristics to the last section. 

Figure 4.1: We show here a model obtained when a 

solar wind velocity of 200 km/sec is used to invert the 
data. Also shown are the standard deviation and resolution* 
* of the parameter corrections when 2, 3, and 4 eigen- 
vectors are used to form the resolution matrix R„. We 


*The elements, of the resolution matrices were rounded to 
increments of 0.2 merely to facilitate their plotting. This 
does not, however, significantly alter their information. 










recall, however, that only 2 to 3 linear combinations of 
parameters were allowed to participate in the iteration 
in order to insure its stability • Perhaps we should not 
be surprised since 5% uncertainties in the data can 
result in a standard deviation of .the parameter correction 
of 60% for four linear combinations of parameters. 

We note that we have the good resolution at shallow 
(< 250 km) , and moderate depths (from about 550 to 750 
km) , but a poorly resolved section between these two 
regions. Moreover, within the linear regime around our 
solution, the data practically are insensitive to the 
conductivity value at depths greater than about 750 km. 

Figure 4.2: We show here the relatively good fit 

to the A . data that occurs when the preceeding con- 
mi n 

ductivity model is used to compute the response. The 
data resolution matrices exhibits a more or less uni- 
form distribution of information among the observations 
with somewhat higher resolution at the high and low 
frequency ends of the spectrum. 

Figures 4.3 and 4.4: The results of a similar cal- 

culation using a solar wind velocity of 400 km/sec, 
are shown here. Again, a relatively good fit to the data 










can be obtained but this higher solar wind velocity results 
in a substantial modification to the shallow conductivity 
structure of the model. The conductivity of the second 
layer was of such a low value that it was practically 
unresolved. A fairly well resolved section at a depth -of 
150 to 350 km was followed by a poorly resolved dip in the 
conductivity structure. Due to the possible importance 
of such a conductivity minimum, we examine its relevance 
in more detail in the next figure. 

Figure 4.5: The middle diagram exhibits two models with 

different layer thicknesses than the preceeding ones, but 
computed with the same solar wind velocity. Their res- 
pective fit to the data is shown in the top diagram. Though 
the model with the more pronounced minimum does give a 
somewhat smaller least squares residual to the data, 
their computed responses differ only at the high frequency 
end of the spectrum. But since the response in the hihg fre- 
quency range is quite sensitive to the source parameters 
we can not consider one model more ade^ 

quate than the other. This is emphasized, by the 
200 km/sec model shown in the bottom diagram which has a 
monotonously increasing conductivity with depth and 
actually gives an exact fit to the A . data values.: 
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7’ 

We recall that in the Inversion we assume that the 
wave normal is parallel to the solar wind velocity. 
However, the average value of the data probably is 
determined from the superposition of a spectrum of fluctu 
ations with significant power in waves with normals at 
substantial angles from the direction of the solar wind 
velocity. We recall from Figure 3.13 that an increase 
in this angle lessens the response at high... frequencies 
and that a similar reduction occurs when response is com- 
puted for decreasing values of the solar wind velocity. 
Thus, even though 400 km/sec is near the observed average 
solar wind velocity of 350 km/sec,, we cannot surmise that 
it is a more adequate value than 200 km/sec to use in 
an inversion which assumes the wave normal to be parallel 
to the solar wind velocity. 

Figures 4.6 and 4.7: We show here three models, 

numbered 1, 2, and 3, with which we attempt to fit 

respectively the A . , A , and A data values obtained 

* mm y z 

at the Apollo 12 site. In all three cases, a solar wind 

velocity of 300 km/sec is used in the inversion. The 

resulting data fit the data shown in Figure 4.7. Notice 

in this last figure that the A value, though generally 

z 

higher than the A -value at most frequencies, is smaller 
than Ay at 0.04 Hz. As we have repeatedly stressed, the 
rollover in the response at high frequency is quite 
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dependent on the parameters of the source . Thus , the 

relatively bad fit to the A data at high frequencies 

is probably symptomatic of our inadequate knowledge of 

the source parameters. Returning to Figure 4.6, we 

note that the conductivity model that attempts to fit 

these values shows a conductivity level reaching 
-2 

nearly 10 mho/m at depths greater than 550 km. How- 
ever, we shall show in the next section that if we toler- 
ate a somewhat worse fit at high frequency, the low 
frequency data can be matched by a model having a sub- 
stantially lesser conductivity at these depths. Since 
we are probably not justified with the data at hand to 
require a close fit to the high frequency response, the 
argument for a high conductivity at great depth becomes 
rather tenuous . 

Before proceeding to our study of conductivity 
structures subject to a given temperature model, let us 
discuss the well resolved features of models 1 and 2. 

We note from Figure 4.6 that at depths between about 50 
to 250 km the conductivity of model 2 is slightly larger 
than that of model 1, This conductivity difference is due 
'to the higher values o'f the A data with respect to the 
A min < ^ ata at high frequencies. However, the approximate 
equality of A m j_ n and A^ at low frequencies requires about 
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the same average conductivity-thickness product for 
both models. Consequently, at moderate depths 
(250 - 450 km) model 2 must exhibit a lower conductivity 
than model 1 to offset the higher conductivity of model 
2 at shallow depths (50 - 250 km) . 

4 . 4 ‘ Inversion of a L-ayered Moon Constrained by a 

Temperature Model 

We would like to know to what extent a given con- 
ductivity profile for the moon will constrain its com- 
position and internal temperature distribution. A 
customary way to study this question is first to assume 
a composition based on independent geophysical and geo- 
chemical evidence. Then, from laboratory measurements 
of the conductivity-temperature relationship of the 
assumed mineral assemblage, deduce the temperature dis- 
tribution from the electrical conductivity model. 

For most types of common reocks and minerals, the 
conductivity-temperature relationship is depicted by a 
series of connected straight lines on a log a versus 
the inverse of the absolute temperature. Each of these 
straight lines can be described conveniently by a relation 
of the type : 

-E '/kT 
= a e o 
o 


a 


(4.4.1) 
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where o and E are called respectively the con- 
o o 

ductivity prefactor and activation energy and where k 
is the Boltzmann’s constant, with the temperature T 
in degrees Kelvin. 

Measurement of relations of the type 4.4.1 often 
are subject to large uncertainties, which arise for various 
reasons, such the fact that minor constituents of the 
mineral assemblage and the oxygen content in the atmos- 
phere under which the conductivity of a rock sample is 
measured can disproportionately affect the bulk con- 
ductivity of the sample. Recent experimental investi- 
gators, especially Duba (1972a), have clarified this 
range of measurement uncertainties and some of their 
probable causes. By measuring the conductivity of 

different olivine samples, with essentially the same 

3+ 

fayalite content, but with differing amounts of Fe 
and under atmosphere with widely different oxygen con- 
tent, Duba was able to show that the temperature estimate 
from the conductivity of the earth's upper mantle might 
be in error by as much as 700 °C. Thus, any attempt to 
deduce an estimate of the temperature inside the moon 
is a hazardous process indeed. Nevertheless, several 
such attempts have been made recently and we listed some 
of them merely to emphasize the range of uncertainties 
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involved. These temperature estimates were aimed especially 

at determining the maximum temperature reached in the 

lunar interior. Since these estimates were based on an 

. . -3 

assumed conductivity level of about 10 mho/m reached 
between 500 and 700 km inside the moon, we would deduce 
essentially the same maximum temperature range from the 
results of the last section. 

Sonett et al. (1972) , using essentially the con- 
ductivity-temperature relationship of olivine from 
England et al . (1968) , have estimated that the maximum 
temperature reached is between 800 to 1000 °C in the 
upper 700 km of the moon. But, based on measurements 
on olivine and pyroxyne in an atmosphere with very low 
oxygen content, a condition believed to occur in the 
lunar interior, Duba et al . (1972b, 1973) were able to 
raise the preceeding estimate to values in the range 
1100 to 15Q0°C. However, Tolland (1974), following a 
suggestion by Ringwood et al . (1970) that pyroxenite 

might be a correct choice for the lunar mantle composition, 
conducted some measurements on a. synthetic lunar pyroxenite 
sample and obtained temperatures in the range 550 to 
680°C. We should point out that Duba et al . (1973) also 
measured the conductivity of a synthetic lunar pyroxenite 
in an atmosphere of very low oxygen content. The con- 
ductivity level of 10 mho/m was reached at about 1000°C 
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but they considered the uncertainties in their measure- 
ments too large to warrant any inference to be based 
on them. 

Some of the data usen in these estimates are shown 
in Figure 4.8. The Fe 3+ free olivine data is from Duba 
et al. (19 72b) and was measured under a pressure of 8kb. 
The synthetic lunar pyroxenite data is from Tolland (1974) 
and the olivine from England et al. (196 8) . Also shown 
are the conductivity-temperature relationships of a 
lunar crystalline rock from Schwerer et al. (1971) 
and two curves labelled H and L that will be discussed 
shortly . 

Due to our present uncertainties in our knowledge 
of the conductivity-temperature relationship of the 
material likely to form the bulk of the moon it might 
be useful to inquire to what extent given temperature 
models for the moon constrain this relationship. One 
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idea behind this question is that the answer would auto- 
matically satisfy the impressive set of constraints that 
can be used in a calculation of the lunar thermal histo- 
ry. (see for example Solomon et al . , 1973). Though these 
cons taints help narrow the range of values that the 
temperature might attain in the lunar interior, there 
is still some controversy about the role of solid state 
convection in the lunar thermal history. This problem 
is still being ..investigated actively 
(Cassen et al . , 1974). But ; due to the uncertainties in 
the viscosity temperature relationship of the material 
inside the moon and the time scale associated with signi- 
ficant heat transfer, by solid state convection, no defi- 
nite conclusions have been reached yet. Never- 
theless, some general conclusions can be reached 

for a rather wide family of temperature model. To get 
a qualitative idea on how this comes about, let us 
consider two temperature models proposed by Tokzoz et al . 
(1972) which were calculated using different 

uranium concentrations ( u m Figure 4.9) . We note that 
in these models that the temperature at depth 200 to 
700 km increases by more than 6Q0°C. Now let us assume 
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that the conductivity at these depths is characterized 
by a conduction mechanism with activation energy of 1 ev, 
an energy which is typical of the range of values for 
the materials shown in Figure 4.8. Then, the electrical 
conductivity would increase by more than 7 orders of 
magnitude between 200 and 700 km and this increase would 
flatly contradict our conductivity estimates of the last 
section . 

Two simple models can explain this discrepancy. 

Either the activation energy is much lower than 1 ev 
or the temperature gradient at these depths is much 
smaller than the one shown in Figure 4.9. Although more 
complex explanations may exist, one of which shall be 
studied briefly later on, we shall concentrate mainly on 
the simpler alternatives since they might be of a more 
general interest. 

In order to find the two parameters in Equation 
(4.4.1) that satisfy both the magnetic data and a given 
temperature model we can use an iteration similar to the 
one of the last section. This time, however, the log- 
arithmic parameter correction must be found through the 
following matrix equation: 

1 -E o /kT 

B ) - E o /k7 2 Alna o 

! AlnE o = P (4.4.2) 
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where B and p are given by Equation (4.2.4) and T. is 
the temperature in the middle of layer i. 

We proceed now to detail the results obtained 
reserving discussion , for the last section. 

We show first, in Table 4.1, the values of the 
parameters obtained for different values for the solar 
wind velocity when the moon is characterized by the low 
temperature model depicted in Figure 4.9. The corresponding 
conductivity models and the fits of their responses to 
the A . data are shown in Figure 4.10. 


Table 

(Data: 


4.1 


A min^ 


V sw ( km / sec ) 


a Q (mho/m) 

E 0 (ev) 

200 


O.lSxlO" 2 

0.093 

300 


0.31xl0” 2 

0.14 

400 


0.47xl0" 2 

o 

• 

CO 


The set of parameters shown in Table 4.2 were ob- 
tained again by using the low temperature model of 
Figure 4.9, but this time an attempt was made to match 
our different sets of data, keeping the solar wind 
velocity at a constant value of 300 km/sec. 
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Table 4.2 


(V = 300 km/sec) 
sw 


Data 

a 

o 


E o 

A . 

0.31X10 -2 


0.14 

mm 

A 

0.17x10 1 


0.092 

y 

V 

0.99xl0“ X 

* 

0.16 


We show on the bottom diagram of Figure 4.11 the 

conductivity models obtained when the and A^ 

data values are used in the inversion. The fit to the 

A and A data is shown in the top diagram, 

y z 

From the last two figures, we note that the computed 
responses do not give a very good fit to the data at the 
high frequency end of the spectrum. But, due to our un- 
certainties in the source parameters, we can hardly rule 
out any of these models on the basis of this misfit. 

Notice also that to fit the A data at low frequencies 

_ 3 

we require only a conductivity of the order of 10 mho/m. 

We 'show in Figure 4.12 the standard deviation and 
resolution of the parameter corrections associated with the 
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Giving Best Fits to A n , A^, A z Data 
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three models of Figure 4.11. If we compare these reso- 
lution matrices with the ones obtained in the preceeding 
section, we note a significant redistribution of the 
resolving power. The new distribution exhibits gradually 
increasing resolution with depth. This redistribution 
varies mainly because the temperature model forces a 
lower conductivity at shallow depth than that required 
for a better fit at the high frequency response. This con- 
straint prevents these conductivities from shielding the' 
lunar interior from induction over a substantial part 
of the frequency range of interest. 

In order to study the dependence of our two parameters, 
°o and E o on different temperature distributions, we con- 
sidered, first the high temperature model shown in Figure 

4.9. The best fit to the A , data, for V = 300 km/sec, 

mm sw 

was obtained for a - 0.18xl0~ 2 mho/m and E =0.12 ev. 

o o 

This pair of parameters pertains to the model labelled H 
in Figure 4.8 whereas the curve labelled L refers to the 
model obtained above using a lower temperature distribu- 
tion but the same A . data and source parameters. The 
bottom scales in Figure 4.8 indicate the depth at which 
a given temperature is reached in each model. Each set of 
crosses indicates the conductivity at depths common to 
both models. Thus, except at fairly shallow depths, where 
there is poor resolution, the conductivity of both models 
is essentially the same . However , somewhat different para- 



meters are required to adjust for the difference in temperature 

#. 

A more comprehensive view of the dependence of our 

two parameters on the temperature distribution can be 

gained by considering the more complete though less 

structured set of models shown in Figure 4.13. In the 

inversions, we have assumed a solar wind velocity of 

300 km/sec and used the A . data. The values of a 

min o 

and E obtained, except for models 1, 2, and 3, are 

o i 

shown in Table 4.3 and they give a fit to the data com- 
parable to the ones above. 



Table 4.3 


Model 

(mho/m) 

E 0 (ev) 

4 

0 . 45x10 "" 2 

• *■ / 

0.24 

5 

0 . 28x10-2 

0.17 

6 

-J*'' 

0 . 2 5x10 ~ 2 

0.13 

7 

0.31xl0“ 2 

0.12 

8 

0 . 2 9x1 0” 2 

0.14 

9 

0.36xl0" 2 

0.13 

10 

0.46xl0~ 2 

0.13 


Models 1, 2, and 3 are exceptions since the use of 

both o and E in our inversions would lead to an unstable 
o o ... 

iteration. Instead we had to use a general inverse opera- 
tion with a value of £ that essentially kept fixed 
in, the iteration and forced the activation - 
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Figure 4.13 





energy to assume the value 


E = kT n In — 
o 1 cr^ 


(4.4.3) 


where is the temperature at depths greater than 2 00 km, 
a Q the initial value of the conductivity prefactor used 
in the iteration, and cr-^ the average value of the con- 
ductivity that must be met at depth greater than 200 km in order 

for our response to fit the data. Clearly then, the only parameter that can 

be inferred from such results is which had a value of 
6.2x10 ^ mho/m for model 2 and was within 5% of this value 
for models 1 and 3. The fit to the data for one of these 

conductivity models is shown in Figure 4.14 and is seen 

to be very good. 

Before closing this" section, let me add a final re- 


mark : 

It is interesting to note that the conductivity obtained 
at relatively shallow depths inside the moon is of the 
same order of magnitude as the one found in the earth's 
crust. However, the conductivity in crustal regions of the 
earth is determined mainly by the water content in the 
rocks. Though the majority of the present geochemical 
and geophysical evidence' from the Apollo mission points to 
a very anhydrous moon, we shall investigate 
vestigate how well the data can be fitted if the moon is 
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Figure 4.14 




assumed saturated with water (except for a non-conducting 
crust of 60 km). To examine this model, we have assumed the 
conductivity inside the moon to be related to the tem- 
perature through a relation of the form 


a = 0 o e 


- Eo /kt 


a e 
w 


0.066/kT 


{4 


The first piece of that expression represents the 
temperature dependence of the material that eventually 
dominates the conductivity at the high temperature pre- 
sumably reached deep inside the moon. The second term 
accounts for the presence of water. The factor 0.066 ev 
in the exponent of this expression was determined as 
follows: 

1. The temperature model inside the moon was as- 
sumed to be the high temperature model shown in Figure 4.9 

2. The pore pressure was assumed to be the litho- 
static pressure and was calculated by assuming a uniform 
moon of density equal to its mean value. 

3. To model the porosity dependence on pressure, we 
used the conductivity measurement made in the laboratory 
on a Westerly granite subjected to pressure up to 10 kb 
at 20°C. 

4. The pressure and temperature dependence of the 
water solution was assumed to follow the laboratory 
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measurement on a 0.01 molar solution of NaCl. 

This type of calculation is very similar to the one 
made by Brace (1971) . who examines the conductivity that 
would be found deep in the earth's crust if it were satura- 
ted with water. The data used in our calculation can be 
found in this paper and the references quoted therein. The net 
result for the moon is a decrease of conductivity with 
depth and is modelled quantitatively by the factor 0.066 
in Equation (4.4.4). 

The Equation (4.4.2) can be easily extended to in- 
clude the inversion of the three parameters o , E and a . 

o o w 

However, in order to insure the stability of the iteration, 
we had to introduce a value of e that would effectively 
permit only two linear combinations of these parameters 
to participate in the solution. Though most of the re- 
solving power was concentrated in E and a , the a para- 
meter had also some, power and did vary in the iteration. 

The model obtained together with its fit to the 
data are shown in Figure 4.15. The bottom scale repre- 
sents the temperature and lithostatic pressure reached at 
a given depth for our assumed temperature model. The best 
fit to the data was obtained for o q ~ 0.054 mho/m r 
Eq ” 0-49 ev and = 0.35 x 10 ^ mho/m. We note that 
the match to the A . data is rather poor for this set of parameters. 
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Figure 4.15 
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In fact, a much better fit could be obtained if water were 
absent altogether and only and E Q were allowed to par- 
ticipate in the iteration. Our inability to obtain a 
good fit from a water saturated moon is indicative that 
the porosity decrease with pressure poorly met the 
variation of conductivity with depth required by the 
measured response. 
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4.5 Conclusion, with Suggestions for Future Work 

The various conductivity models obtained in the last 
two sections. are representative of our present 

uncertainties arising from our lack of knowledge of the values of 
various parameters that influence the data. Yet all these 
models seem to possess some common characteristics. 

We note first that a conductivity level of 10 ^ mho/m 

seems to be consistently reached in the upper 150 to 200 km 

of the moon. Moreover, the conduction level does not generally 

increase by more than about one order of magnitude between 

200 and 700 kilometers. Also, the conductivity below 

about 8 00 km is generally unresolved by the data. Though our studies 

of the A z data from the Apollo 12 site does not rule out' 

-2 

models in which a conduction level slightly less than 10 mho/m 
is reached in the top 700 km of the moon, vje feel that a lower con- 
duction level is more likely since the A component 

4-r 

lies near the direction of the tangential component of the 
remanent field at this site. Moreover, the low frequency 
data at the Apollo 15 site tends to aqree with the A . 
and A y values from the Apollo 12 LSM and these data yield 

_ "5 

a conduction level consistently less than 2 x 10 mho/m 
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in the upper 700 km of the moon, 

v 

If we attempt to extract the maximum temperature 

♦ 

reached in the moon from the conduction level, we face 
such large uncertainties that we cannot rule out a tem- 
perature either substantially higher or lower than 1000 °C. 

To narrow the allowable temperature range, a great deal more work must 
.be done to obtain and 

understand the conductivity-temperature relationship of 
mineral assemblages likely to form the bulk lunar composition* 


Yet, the fact that the conduction level does not 
vary much between 200 and 700 km inside the moon seems to 
put a rather stringent constraint on either the conductivity- 
temperature relationship or the temperature structure at 
these depths. We have examined this question in the last 
section and found that if a sharp thermal gradient is 
maintained at these depths, a very low activation energy 
conduction mechanism is required. In fact, even if 
the sharp thermal gradient is maintained only in the top 
300 km, an activation energy of only about 0.24 ev 
is required. For models more in accord with the 
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one proposed by Toks<5z et al . (1972) and in which the sharp 

thermal gradient extends to depths greater than 500 km, 
the activation energy required is generally smaller than 
0.2 ev. A low of 0.09 ev can be reached depending 
upon the set of data and source parameters used in the 
inversion. Moreover, these results do not seem very sen- 
sitive to the maximum temperature reached in the lunar in- 
terior. 

Low activation energy conduction mechanisms are often 
observed in terrestrial and lunar* rocks (see, for example 
the low temperature end of the lunar sample a versus inverse T 
curve shown in Figure 4.8). For most cormon rocks, such conduction 
mechanisms are observed only at fairly lew temperatures. Moreover, they 
have an associated conductivity prefactor generally much 

-3 

smaller than the level of a few times 10 mho/m that 
would be required for the rocks inside the moon. Thus, 
there are definite problems with our inferred model parameters 

We note that low activation energy is often 

associated with impurities in the major mineral component 
of the rocks. Another possible interpretation is that 

Preliminary measurements of the conductivity-temperature 
relationship in lunar rocks by Schwerer et al . ( 1971 ) were pla 

gued by thermal hysteresis that was subsequently attributed to a carbon 
contamination of the sample. The implication of these uncertainties 
is that, at the present stage, very few results can be relied on. 
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the rocks contain a small amount of a mineral whose con- 
ductivity is high with a weak dependence on the tern- 
perature. A fairly common terrestrial mineral that pos- 
sesses these • characteristics is magnetite. Above 119 °K, 
its conductivity is of the metallic type and has a value 
of about 2 x 10 ^ mho/m (Miles et al . , 1957) . The con- 
ductivity of hematite doped with a small amount of titanium can also 
exhibit a weak dependence on temperature and a fairly high conduction 
level (Morin , 1951) . 

The fact that our low values of E and moderate value of a 

o o 

have not been observed in the Apollo samples might be the result 
of a significant differentiation that has occurred during 
the moon's history. This differentiation could have 
depleted the lunar surface rocks of a highly conducting 
yet minor mineral component. The enrichment of this com- 
ponent at moderate depth could account for our inferred 
characteristics . 

If a minor mineral dominated the conductivity at 
moderate depth inside the moon we cannot determine 
from the magnetic data the major mineral component 
at these depths. We might even have difficulty predicting the con 
centra tion of the minor mineral since the bulk. 

conductivity probably depends on the connectivity of its 
conduction paths , a factor over which we have little 
control. 

i pRTCTN A B PAGE IS 
m POOR $0ALKrei 



The dif f erentiation ' migfit partly explain a dif- 

ficulty that could result if the near surface rocks have 
the same conductivity parameters that we inferred from the LSM data 
at moderate depth. From Figure 4.8 we note that the conductivity that we 

would extrapolate for rocks near the surface would be of 

-5 

the order of 10 mho/m. But, if such a high value were realized, 
a strong toroidal H magnetic field would be expected, 
contradicting our observations. However, differentiation 
might not be required to explain the absence of this field. 
All we need is that the moon be covered by a layer of 
resistivity-thickness product greater than about 
5 x 10^ ohm-m~km. On the basis of high frequency electro- 
magnetic data (136 Mhz to 75 Ghz) , Strangway has inferred 

the existence of' a powdered rock layer of conductivity in 
— 1 3 — 1 6 

the range 10 to 10 mho/m near the lunar surface. 
Therefore, the thickness of this dust layer need only to be 
less than about 1 meter to explain the absence of perturba- 
ion. Also the intense scattering of seismic waves that is 
observed on the moon has been widely interpreted to imply 
the existence of poorly consolidated material near the 
surface. Evidently, if the moon is dry, the increased 
porosity will break up the conduction path and might result 
in a greatly reduced conductivity, even if the rocks near 
the surface had the same composition as those at moderate 
depth . 
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An alternative explanation for the absence of strong 

* 

variation in the electrical conductivity at depths below 
200 km would.be to assume that there is no significant 
thermal gradient at these depths. We have examined this 
possibility in the last section and we found that indeed 
a reasonably good fit to the data can be obtained in that 
case. Such thermal models do not pose any significant 
constraint on the activation energy of the material at 
these depths. Kuckes (1974) has also recently investigated 
this question. He assumed at the onset that the activation 
energy is higher than about 1 ev and searched -for 
temperature model yielding a good fit to the magnetic data. 

He also found that only the top 200 km can maintain a strong 
thermal gradient in order to meet his constraint on the acti- 
vation energy. 

Temperature models with small thermal gradients below 
200 km have been proposed by Tozer (1972) . From considera- 
tion of material rheology under physical conditions likely 
to be found inside the moon, he concluded that a convecting 
core probably exists within 200-300 km of the surface. 

Bypass uming the Nabarro -Herring creep 

to be the dominant deformation mechanism and by inferring 

21 

a viscosity of 10 poise in the convection region, he was 
able to establish a lower bound of 5. x 10 mho/m on the 

electrical conductivity that would be permissible inside 







- 160 - 


the moon. Since we deduce a conductivity level of about 
6.2xl0 - ^ mho/m when such a temperature model is assumed,' 
we cannot reject his hypothesis on the basis of the mag- 
netic data alone. However, the lower bound on 
the conductivity quoted by Tozer is not 
very well established in view of the fact 

that the Nabarro-Herring creep is only one member of a set 

of creep processes that could have participated in such a 

21 .. 

convection. Also, a viscosity of about 10 poise is 
often associated with regions of the earth's upper mantle 
where creep is presumed to be an important deformation mechanism 
(see, for example, Weertman, 1970). But these regions of 
the upper mantle are characterized by values of the elec- 
trical conductivity several orders of magnitude greater 
than 'that we would infer for the moon. However, the different 
physical conditions in the earth from those in the moon 
might perhaps account for some of these discrepancies . 
Nevertheless, at this point, a widely recognized objection 
to the calculation of Tozer is the implication of a rapid 
cooling of the lunar interior after the moon's formation. 
This cooling rate would seem to contradict the substantial 
period of time ( 0.7 to 1.4 by) during which the basalt 
that flooded the maria was kept near its melting temperature 
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in the early history of the moon (Ringwood et al. , 1970) . 
Until this objection can be countered, the results of 
Tozer will remain inconclusive 

Thus, our results seem 
to lead to somewhat puzzling conclusions. 

A substantial effort should be devotea in the future to 
either confirm or deny their validity. We have pointed 
out in the text various directions that can be followed in 
order to assess the accuracy of the theory and the possible 
noise contribution to the data. 

Let me just add one more suggest- 
ion. Future analysis should include in the inversion 
an accurate value for the position of the LSM, che solar 
wind velocity and the distribution of the wave normal to 
the incident fluctuations. The value of the solar wind 
velocity probably can be obtained from the various instruments aboard 
spacecrafts in the solar wind or from the solar wind spec- 
trometer at the Apollo 15 site. The small remanent field 
at this site probably does not influence the velocity. The 
distribution of the wave normal may be more difficult to 
reconstruct. However, the situation is not completely hope- 
less in view of the preliminary results of Daily (1973) which 
tend to indicate that during a good Alfven wave regime, the 

ORDINAL 

0F fOOE QUALTIm 
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wave normal tends to align along the ambient magnetic field 
direction. If this is the case, then, during such periods 
the problem would reduce to that of finding the average 
direction of the magnetic field over a period appropriate 
for the frequency range of interest. Alternatively, we 
might be able to extract some of this information from 
the strong dependence of the high frequency response in the 
downstream cavity on the value of these parameters. We 
visualize, for example, a situation where Explorer 35 and 
the LSM are on the upstream side of the moon when a 
subsatellite is in the downstream cavity. The high fre- 
quency response at the subsatellite possibly could be used 
to extract some of the required information which would 
in turn be fed to the inversion of the LSM data. However, 
it might be impractical to extract useful information 
from the subsatellite due to its relatively short orbital 
period and consequent short sampling time, coupled with our 
inferred strong dependence of the high frequency signal in 
the void on the distance from the plasma-vacuum interface. 
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Appendix I 

The Noise from the Solar Wind Dynamic Pressure Fluc- 
tuation at the Apollo 12 Site 


We can not estimate accurately the effect of the 
solar wind dynamic pressure fluctuation on the remanent 
field around the Apollo 12 site due to our incomplete 
knowledge of the structure of the remanent field it- 
self. Moreover, the solar wind data necessary to back 
up such an estimate are generally not available. How- 
ever, we may be able to extract the order of magnitude 
of the solar wind pressure effect from some empirical 
relations* that have appeared recently in the literature. 

We note first that Dyal et al . (1972) induced an 

empirical relation between the change in the magnetic 
field at the site and the solar wind dynamic pressure. 

They have expressed this relationship as follows: 


* 

The cgs and Gaussian system of units is used in this 


note . 
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(AP-I.l) 


and B is the total field at the site. B_ is the field 

E 

k 

measured by Explorer 35, B is the unperturbed remanent 

field and p is the solar wind density. Hourly averages 

of 5 and B_ were used to estimate AB so little con- 
E 

tamination from the induced field is expected. 

At first glance Bguation (AP-I.l) appears to be an outgrowth 
of the conservation equations. However, we can make such an inter- 
pretation only if. the 

square of difference of |S - S^| and [§ s | had been 
related to the dynamic pressure (see, for example, 

Siscoe et al. , 1971). 

During a period when the solar wind magnetic field fluctuations 
were predominantly Alfven wave, Belcher et al. (1971) found the fol- 
lowing relation between the solar wind magnetic fluctuation, 6B , and 
* > 

the fluctuation in the solar wind velocity, 5V 

2 sw: 

aS sw = ±4.6 p 1/2 6^ sw (AP-I.2) 

In order to relate these two empirical equations, let us 
assume that the modulation in the dynamic pressure due to 

■u 

velocity fluctuations along the average direction of the 
solar wind velocity is associated with a perturbation 


A§ • AB 


0 = 0.01 pV -V , 

8tt sw sw 


where 

. v 


*► ->• *>■ 
AB = B - B E - B s 
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6B , of the remanent field. If we linearize 
mod. 

Equation (AP-I.l), that is, we assume AB = AB Q + 

— y ~y Q -'V t 

and V = V . + 6V . and if we collect the first order 

sw sw sw 

terns we obtain 




= 0.01 pV • <$v 

sw sw 


If we define the direction of the unit vectors a^ and 
a v parallel respectively to AB q and V^ w and using the 
equation obtained by collecting zero-order terms, we obtain 


L b ^ B mod 


n _ 1/2 

^ 0 . 5 p a *oV 

V sw 


which, using Equation (AP-I.2), becomes 


0.09 a • 5B 

v sw 


Consequently we expect a ratio of the maximum to minimum 
tangential amplification factor smaller than about 1.1. 
However, Sonett e t al . (1972) have observed a ratio equal 

to about 1.5 at 0.005 Hz. It would seem difficult to 
interpret such a high value on the basis of our rough 
amplitude estimate. There are two observations that also 
seem inconsistent with an interpre tation in terms of 




pressure fluctuation noise. One: Sonett et al . found 


that the direction of maxima in the amplification factor 
does not align with the direction of the tangential mag- 
netic field component for frequencies below about 0.002 Hz. 
Second.: they also found two directions of maxima in 

the amplification factor at 0.04 Hz. However, more 
experimental and theoretical work is needed before we can 
properly assess the contribution of pressure fluctuations 
in the anisotropic response observed at the Apollo 12 
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Appendix II 

The Field Inside the Moon 


The solution of Maxwell's equations inside a layered 
sphere has been discussed frequently in the literature. 

We shall use here the results and notation of Stratton 
(1941) and shall content ourselves to describe the quan- 
tities L and T that permit us to write the field as in 
n n 

Equations (2.3.1) . 

If we assume the moon to be formed of s_ layers 
bounded by the radii r^, r • • • r r where r g = R^, then 
and can be expressed as follows : 


L n = 1 + D 2 /D 1 (Ap. II. 1) 

where 


D 1 

= M (k r )M ^(k r , ) M (k , r -,)... 

j n (k l r l> 

D 2 

n s s n s s-1 n s-1 s-1 

k l r in (k l r l ) 


and where the 2x2 matrix M (z) is given by 


M n (z) 


j n (z) n n (z) 
zj * n (z) zn^(z) * 


Similarly we have 


T 

n 


V F 


2 


where 


(Ap. II. 2) 



■-1 


N n (k s r s > N n (k s r s-l> K n (k s-i r s-l ) 


k l r l3n (k l r l> 

K ij-tk 1 r 1 j n (ki r 1 )] 1 
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where the 2x2 matrix N n (z) is given by 

z j ( z ) zn {z ) 

. , , * J n n 

N n (z> = 

--[zj (z))' — [z.n (z) ] ' 

z J n z n 

v 

1/2 

where = (icoy 0 ^^) is the propagation constant 
of the layer 

and n n (z) and j n (z) are the spherical Bessel functions and 
the prime stands for the derivative with respect to the " 
argument. The Inverses of M n and can be found easily 
and simplified by using the Wronskian re- 
lationship for spherical Bessel functions. 
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APPENDIX II 

Response in the Void Cavity to a Magnetic Discontinuity 
in the Solar Wind 


We would like to calculate the field in the cylindrical 
void region when a magnetic discontinuity is convected 
in the polar wind and compare the result with the 
measurements exhibited in Figure 2.4. During these 
measurements, Explorer 35 was in the void region behind 
the moon and both the moon and Explorer 33 were in the 
free streaming solar wind. We shall assume that Ex- 
plorer 35 was sufficiently far downstream that its mag- 
netic field measurements were uninfluenced by the con- 
ductive lunar interior [see Sonett et al . (1970) , for 

a discussion of this assumption] . We shall also assume 
the incident solar wind magnetic field to be homogeneous 
in planes subtending an angle tf; with the direction of 
the solar wind velocity. However, we have only one 
relation to determine ip and V sw , namely, the delay in ar- 
rival time of the discontinuity at both satellites 
(- 15 minutes) . From the difference in the X and Y 


coordinates of each satellite, (AX = 22R , AY = 48R ) , 

Hi ill 

-we obtain 

AT = 900 sec = tanj, (AP.II 

v sw 

If we assume ti) = 0, we obtain V - 150 km/sec. This value is 

sw 

substantially lower than that of the solar wind velocity (=350 km/sec) . 
Nevertheless, first we make a calculation assuming this low solar wind 

velocity and then present a brief heuristic argument to show that even 

for very high solar wind velocities, the predicted dilatation in the 



rise time of the Z exponent ranains of the same order of magnitude 
as the one observed. 

If Of then, in the XZ plane (see note, p. 53 and the table 

of Figure 2.4) field predicted by Equation (2.4.1) is 

ik , ,z 


H„ = H = 
Z p 


H ox (M)I i (k // p)e 

ii< k // a > 


// 


and 




H„ = 


H * = 0 


ik 


-H = 
z 


- H ox (w)I l (k //' p)e 




"iV 


(AP. II. 2) 


where the Fourier transform of H (w) is assumed to be the 

ox 

signal observed at Explorer 33. The Fourier transform of 
these expressions have been evaluated and are 
presented graphically in Figure AP-III. The 

dilatation in the rise time of the Z is about the same as 
the one observed and the peak in the X is also along the 
same direction as observed. 

Now, let us consider the case of a high solar wind 
velocity. From Equation (AP.I 


AT 

AY 


tantp 

V 

sw 


If the solar wind is allowed to reach very high values, but 
the above ratio is kept constant, we obtain from 
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Equation (2.4.1) 


ik /, 2 J (aw) 

H„ = H « H (w) e 7/ l ~ 

Z P ox' ; . aw 

rn ± ^ 3 / • * • 


where 


ATR 


a - 


M 


AY 


We can easily show that the min contribution to the 
dilatation in the rise time of the Z component will come from the term 
m = 1. This term has a Fourier transform given by 

~ [1 “ (t/a) 2 ] 1 ^ 2 for |t| 1 a 


t A a 


The convolution of this Fourier transform with the source field 
adds about 2a - 10 sec to the rise time of the Z component, thus yield- 
ing a total value for the rise time equal to about a third of what is 
observed . 

Since actual values for ^ and V for the discontinuity are 
probably between these two limits* described above, we conclude that 
the surface wave solution of the field in the void does predict, at 
least qualitatively at this time, the observed dilatation in the rise 
tirra of tile Z component and the small peak in the X component . 


4* i.. 

Negative value of can render V smaller than 150 km/sec, 

sw 

but such a case must be considered improbable. 
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Appendix IV 

The Multipole-cylindrical Waveguide Mode Representation 
of the End Effect Field 


Let us consider the free space magnetic field of 
a magnetic multipole (neglecting displacement current) . 

This field can be written as follows: 

.n+2 

m c i t-i -v 

(Ap. IV. 1) 


H = V x V x [t° 

M nm n 


P™(cos0) { Sln } m<j)a ] 
n cos r 


In order to satisfy the condition that the magnetic field 
has a vanishing normal component at the boundary of the 
cylinder, we must find a series of spherical TF 
modes, regular at the origin such that the sum of the 

normal component of the magnetic field of the TE modes and that 
of the irtultipole add up to zero at the boundary. We can ex “ 
press this field representation as follows: 


->■ 

H. 


I V x V x (a 
p=m 


p+i 

r r rsin, . ■> 

„ _ — P COS0) { } mc|)a 

nmp „p-l P cos f r 


(Ap. IV. 2) 


The coefficients in AP. IV. 2 can be 
found by using the following dual pair of identities: 

co 

q 1 ' K m (tp)cos (tzjdt = £.(£~ m ) i r " ^”^P™(cos8 ) (Ap. IV. 3) 


if m + a is even 
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/°° t^K (tp j sin (tz)dt = ( cose ) 

0 m 2 x. 


if m + ! is odd 


00 „ / . A Z 

. _ /4 _ > „ . Jl-m (tr) m, Q x 

cos tz V tp) = , £ 1 11+ioT p Jl (cos0) 

J6=m 


if £ + m is even 


(Ap. IV. 3) 


sin tz I (tp) = Z 

&=m+l 


, 8,-m-l (tr) 2, m, n . 

1 u+m) r P * (cos0) 


if £ + m is odd 


where I and K are the hyperbolic Bessel functions, 
mm 

These identities are consistent with the Morse et al . 
definition (p. 1325) of the Legendre functions and are 
derived by Cooke (1956, 1962). 

Let us assume n + m is even in Equation (Ap. IV. 1) 
and consider the following cylindrical TE magnetic 
field: 

H x = V x V x iff a z (Ap. IV. 4) 

CO 

ib = b a f t n (tp) sin(tz)dt 
Y nm q m r 

a n+2 

, a „_2_ n nm m --sin-, , 

D nm ” it .n-m+2 (n-m) ! l cos 


where 
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By using the first pair of identities in 
Equation (Ap. IV. 3) and making liberal use of Legendre 
function identities (in addition to applying the operator 
VxVx to cancel some terms not contributing to the field) , 
we can show that =-fi . 

In order to cancel the normal component of the 
magnetic field at the boundary of the cylinder, we can 
add to r the following TE magnetic field: 

§2 = V x V x $ a z (Ap.IV.5) 

where 

* u o ,Tn-l K m (tR M ) rL _ % ^ 

* = -b nm f Q t I m (tp)sm(tz)dt 


By following a process similar to the one described 
above but this time using the second pair of identities 
in Equation (Ap.IV.3) and integrating term by term, we 
can show that = H c , provided that we set 


nm 


-y e y e h N (n+ 0 
nm mp nmp m r 


(Ap. IV. 6) 


where 


.a . p-m 
t n i. 
nm 


nmp 


_____ /S+Ri i . 

(n-m) ! (p+m) i (n+ptl) (p+1) s 2 1 { 2 


The notation for is explained in Equation (2.4.1) 

and the constants N (21) are defined in Appendix V. 

The spherical TE electric field associated with 


H m +H c does. not have a vanishing tangential electric field 
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component on the boundary of the cylinder. In order to 
meet this boundary condition, some spherical TM modes 
must be added. They are represented by the po- 

tential in Equation (2.4.3), whereas the electric field 
associated with the electric multipole is represented by 
the potentials in that equation. Since the tech- 

nique used to find these potentials is very similar to 
the one described above, the details of the derivation 
are omitted. 

An alternative representation of the field in terms 

V 

of cylindrical waveguide modes can also be obtained. We 
note that the sum of the two cylindrical TE potentials can 


be given as follows: 


...a ... b nm ” 

vp =^ + $ = _/ t rr/ r j 

nm r 2 I to 


*sintz dt (Ap.IV.7) 

where we have used the definition of hyperbolic Bessel 
functions (Watson, 1966) to extend the integration 
from'- 00 to The integrand is analytic and one-value 

and has simple poles at ±i£ J^(B m ^) = 0. Using 

the residue theorem and noting that the integral vanishes 
on the semi-circle at infinity, we obtain, for z > 0 
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¥ 


CT 

run 


^ a 2 
2nt 

nm>I v . n-2 
(n-m) ! 


^ ^p/R^) e 


“ e m£ Z / R M 




r sin 
l cos 


} m<j> 


(Ap.lV.8) 


The Wronskian relationship for Bessel functions was 
used to obtain this last expression. The electric field 
can also be treated in a similar way and the total field 
is expressed in Equation (2.4.4). 
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Appendix V . 

The Constants M (2p) and N (2p) 
m m ^ 


Several problems in potential theory, involving 
spheres inside cylinders, required the same set of 
definite integrals that we encountered in Chapter II 
(see, for example, Cooke, 1962, and Smy the , 1963, 1964, 

196 8) . Due to the general usefulness of these integrals and since 

only a partial listing appears in the literature, we 

include here a more comprehensive table of their values. 

We define 


Mln (2p) 


2p 


» K„(t) 


m 


Trpip 


1 n t 


0 I (t) 
m 


- 2 (2p+l) r 2p 

- .P!P! £ l m <t) 


dt 


and 


N m (2 P ) 


2 “ t 2 P +1 1 m 2 

— ~ r f — 7T ( =r + 5™ )dt 

tt p.p. o (i * ) 2 t t 3 


m 


= -2 (2p+l) ,°° 11 PK m (t)dt 

Tipipi n 






The above identities can be derived by integrating 
the left hand side and then using the Wronskian re- 
lationship. To evaluate these integrals, we have 
proceeded as follows: 

For p - 30, an analytical expression was used. 
Following a technique devised by Watson (1930), we can 
derive 'the following identities: 
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f 

0 I* 


2p 


dt ■= a S 


(no + 0.5a) 


2p 


m 


n=0 I (na + 0.5a) 


and 


2p-l 

oo f* - 

i „ / n \ v in A 

+ 7T(-1) l —x 


J m+l (o W ) 




cosh ( 7Ta m ^/ G ) 


(2p+l) e 


~™mJ l /0 


cosh (ioP^To) 


/ 

0 


2p 


(I ' ) 
in 


dt = a 


(na + 0.5a) 


2p 


n=0 [I^(na + 0.5a)] 


- TT (-1 ) 


p-m 


oo 

z 




2p+3 


3 m^ /a 


* =1 <4- m2)2j m ( ^£’ cosh 2 (n3 mJ j/a 


- (2p + 


where a n and B „ are the roots of J and J’ respectively 
mi m£ m m 

and a is an arbitrary positive number. 

We used the Royal Society Mathematical Tables (1960) 
to find the roots of the Bessel functions and their as- 
sociated values. The above series were evaluated 
numerically for different values of a. The results 
agree to more than 8 digits. For 50 - p ^ 25, we used 
the asymptotic formula for the product of hyperbolic 
Bessel functions (see, for example, Olver , 1965). In 
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*r 

the overlapping domain 30 - p 1 25, the results were also 
found to agree to more than 8 digits. 

Note: N Q (2p) = M-^ (2p) fJh so the values of 

N (2p) can be found from (2p) 
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Table AP-V : Values of M (2p) 

m 


m 

0 

1 

2 

3 

P 

0 

0.87069013 




1 

1.2354675 

4.7809450 

— 

_ 

2 

1.6471421 

3.0011759 

23.091507 

_ 

3 

2.0014447 

2.9004477 

9.1317676 

105.75794 

4 

2.3046888 

3.0068912 

6.7887287 

30.338664 

5 

2.5708931 

3.1638910 

5.9624110 

18.231357 

6 

2.8105215 

3.3334854 

5.6044969 

13.777881 

7 

3.0303496 

3.5038177 

5.4447761 

11.587752 

8 

3.2347572 

3.6708310 

5.3845295 

10.336636 

9 

3.4266969 

3.8331760 

5.3798713 

9.5556150 

10 

3.6082622 

3.9905271 

5.4 086175 

9.0401532 

11 

3.7810044 

4.1429611 

5.4584810 

8.6879044 

12 

3.9461150 

4.2907110 

5.5221743 

8.4423901 

13 

4.1045347 

4.4340622 

5.5951411 

8.2700663 

14 

4.2570235 

4.5733082 

5.6744117 

8.1498362 

15 

4.4042059 

4.7087308 

5.7579850 

8.0678099 

16 

4.5466033 

4.8405929 

5.8444761 

8.0144936 

17 

4.6846564 

4.9691357 

5.9329054 

7.9831931 

18 

4.8187416 

5.0945793 

6.0225677 

7.9690624 

19 

4.9491838 

5.2171237 

6.1129481 

7.9685147 

20 

5.0762658 

5.3369505 

6.2036664 

7.9788426 

21 

5.2002352 

5.4542245 

6. 2944395 

7.9979687 

22 

5.3213107 

5.5690959 

6.3850554 

8.0242747 

23 

5.4396865 

5.6817015 

6.4753548 

8.0564838 

24 

5.5555358 

5.7921662 

6.5652179 

8.0935769 

25 

5.6690141 

5.9006043 

6 .6545550 

8.1347322 

26 

5.7802615 

6.0071206 

6.7432992 

8.1792806 

27 

5.8894049 

6.1118117 

6.8314013 

8.2266732 

28 

5.9965595 

6.2147663 

6.9188258 

8.2764559 

29 

6.1018303 

6 . 3160666 

7.0055477 

8.3282508 

30 

, 6.2053134 

6 . 4157884 

7.0915504 

8.3817414 

31 

6.3070970 

6.5140023 

7.1768239 

8.4366609 

32 

6.4072623 

6.6107736 

7.2613632 

8.4927835 

33 

6.5058843 

6.7061633 

7.3451674 

8.5499172 

34 

6.6030321 

6.8 002282 

7.4282388 

8.6078979 

35 

6.6987700 

6.8930214 

7.5105320 

8.6665848 

36 

6.7931578 

6.9845926 

7.5922037 

8.7258569 

37 

6 . 8862510 

7.0749883 

7.6731119 

8.7856096 

38 

6.9781015 

7.1642523 

7.75331S0 

8.8457524 

39 

7.0687578 

7.2524256 

7.8328260 

8.9062067 

40 

7.1582654 

7.3395470 

7.9116527 

8.9669041 

41 

7.2466669 

7.4256530 

7.9898074 

9 .0277850 

42 

7 . 3340024 

7.5107780 

8.0673015 

9.0887973 

43 

7 . 4203094 

7.5949545 

8.1441467 

9.1498953 

44 

7.5056235 

7.6782133 

8.2203549 

9 .2110392 

45 

7.5899783 

7.7605837 

8.2959379 

9 .2721939 

46 

7.6734054 

7.8420934 

8.3709073 

9.3333286 

47 

7.7559348 

7.9227685 

8.4452751 

9.3944162 

48 

7.8375948 

8 .0026342 

8.5190526 

9.4554330 

49 

7.9184124 

8.0817141 

8.5922513 

9.5163520 

50 

7.9984131 

8.1600309 

8.6648824 

9.5771729 
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Table AP-V : Values of M (2p) (continued) 


m 

4 

5 

6 

7 

P 





4 

470.30450 

- 

— 

■ — 

5 

105.56602 

2052.4525 

- 

. — 

6 

53.121319 

377.86834 

8840.6117 

- 

7 , 

35.112697 

163.37305 

1378.4996 

37714.041 

8 

26.603247 

9 5.801312 .» 

522.12060 

5097.6836 

9 

21.828804 

65.888357 

274.86722 

1717.1999 

10 

18.847131 

49.917485 

172.85264 

819.69548 

11 

16.845010 

40.311179 

121.43253 

474.42728 

12 

15.429216 

34.038063 

91.938586 

310.38676 

13 

14.388950 

29 .689158 

73.423098 

220.92138 

14 

13.602071 

26.534852 

61.001654 

167.13928 

15 

12.993382 

24.165146 

52.235259 

132.38598 

16 

12.514267 

22.334390 

45.797192 

108.64531 

17 

12.132019 

20.887607 

40.915497 

91.698935 

18 

11.823883 

19.722789 

37.115907 

79.16512 6 

19 

11.573577 

18.770360 

34.093690 

69 .620176 

20 

11.369162 

17.981421 

31.645544 

62.1719 84 ... 

21 

11.201701 

. 17.320699 

,29.631438 

56.238971 

22 

11.064381 

16 . 762177 

27.952234 

51.428997 

23 

10.951919 

16.286283 

26.536051 

47.469861 

24 

10.860167 

15 . 878035 

25.329666 

44.167782 

25 

10.785818 

15.525790 

24.292944 

41.381717 

26 

10.726211 

15.220373 

23.395111 . 

39,006973 

27 

10.679181 

14.954462 

22.612229 

36.964477 

28 

10.642950 

14.722154 

21.925422 

35.193561 

29 

10.61604 4 

14.518634 

21.319634 

33.647023 i 

30 

10.597232 

14.339943 

20.782721 

32.287668 '■ 

31 

10.585482 

14.182799 

20.304792 

31.085845 

32 

10.579916 

14.044455 

19.877717 

30.017667 

33 

10.579791 

13.922602 

19.494759 

29.063704 

34 

10.584467 

13.815280 

19.150289 

28.208011 

35 

10.593396 

13.720820 

18.839577 

27.437396 

36 

10.606103 

13.637789 

18.558615 

26.740860 

37 

10.622174 

13.564954 

18.303990 

26.109171 

38 

10.641250 

13.501245 

18.072780 

25.534527 

39 

10.663015 

13.445731 

17.862467 

25.010296 

40 

10.687192 

13.397599 

17.670873 

24.530809 

41 

10.713537 

, 13.356135 

17.496104 

24.091191 

42 

10.741832 

13.320711 

17.336505 

23.687235 

43 

10.771887 

13.290770 

17.190627 

23.315287 

44 

10.803531 

13.265820 

17.057193 

22.972165 

45 

10.836609 

13.245421 

16.935076 

22.655083 

46 

10.870987 

13.229183 

16.823276 

22.361592 

47 

10.906541 

13.216754 

16.720904 

22.089534 

48 

10 .9 43162 

13.207819 

16.627167 

21.837000 

49 

10 .980750 

13.202097 

16.541356 

21.602294 

50 

11.019216 

. 13.199331 

16.462834 

21,383906 
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Table AP-V : Values of N (2p) 

m 


m 

1 

2 

3 

5 

P 

1 

5.9937671 




2 

4.0668787 

25.388848 

- — 


3 

3.8036566 

10.733735 

111.53052 

— 

4 

3.7744891 

8.0327134 

33.608518 

486.69676 

5 

3.8282347 

6.9789396 

20.480469 

113.35747 

6 

3.9205355 

6,4 638932 

15.477384 

57.879120 

7 

4.0322344 

6.1903803 

12.947072 

38.411860 

8 

4.1536961 

6.0448087 

11.467020 

29.074994 

9 

4.2797371 

5.9742794 

10.522919 

23.780132 

10 

4.4074630 

5.9508685 

9.8862942 

20.447047 

11 

4.5352189 

5.9585102 

9 . 4410020 

18.194926 

12 

4.6620474 

5.9873603 

9.1220915 

16.593983 

13 

4.7873962 

6.0310883 

8.8905584 

15.412127 

14 

4.9109539 

6.0854647 

8.7216496 

14.514087 

15 

5.0325563 

6.1475737 

8.5989470 

13.816190 

16 

5.1521295 

6.2153529 

8.5111607 

13.264142 

17 

5.2696562 

6.2873114 

8.4502900 

12.821307 

18 

5.3851541 

6.3623516 

8.4105192 

12.462138 

19 

5.4986628 

6.4396527 

8.3875273 

12.168323 

20 

5.6102348 

6.5185930 

8.3780431 

11.926418 

21 

5.7199300 

6.5986959 

8.3795484 

11.726346 

22 

5.8278120 

6.6795928 

8.3900771 

11.560418 

23 

5.9339457 

6.7609963 

8.4080735 

11.422671 

24 

6.0383956 

6.8426808 

8.4322930 

11.308416 

25 

6.1412253 

6.9244685 

8.4617292 

11.213919 

26 

6.2424962 

7.0062182 

8.4955611 

11.136170 

27 

6.3422674 

7.0878182 

8.5331131 

11.072720 

28 

6.4405958 

7.1691793 

8.5738251 

11.021558 

29 

6.5375355 

7.2502310 

8.6172294 

10.981017 

30 

6.6331381 

7.3309170 

8.6629333 

10.949705 

31 

6.7274528 

7.4111930 

8.7106049 

10.926452 

32 

6.8205261 

7.4910241 

8.7599623 

10.9102 66 

33 

6.9124023 

7.5703833 

8.8107648 

10.900305 

34 

7.0031236 

7.6492496 

8.8628062 

10.895844 

35 

7.0927298 

7.7276075 

8.9159090 

10 . 896263 

36 

7.1812587 

7.8054455 

8.9699197 

10.901022 

37 

7.2687464 

7.8827556 

9.0247052 

10.909655 

3 8. 

7.3552271 

7.9595327 

9.0301439 

10.921753 

39 

7.4407331 

8.0357741 

9.1361528 

10,936957 

40 

7.5252955 

8.1114789 

9.1926256 

10.954955 

41 

7.6089435 

8.1866478 

9.2494908 

10.975468 

42 

7.6917052 

8.2612829 

9.3066804 

10 .998251 

43 

7.7736073 

8.3353874 

9.3641345 

11.023086 

44 

7.8546750 

8.0489653 

9 .4 2180 00 

11,049777 

45 

7.9349327 

8.4820213 

9.4796304 

11.078153 

46 

8.0144035 

8.5545607 

9.5375841 

11.108057 

47 

8.0931094 

8.6265892 

9.5956244 

11.139351 

48 

8.1710717 

8.6981129 

9.6537189 

11.171908 

49 

8.2483104 

8.7691383 

9 .7118336 

11.205616 

50 

8.3248449 

8 . 835 6 720 | 

9.7699578 

11.240372 
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Table AP-V : Values of N (2p) (continued) 

m 


m 

5 

4 6 

7 

8 

P 





5 

2102.3544 

- 

- 

- 

6 

398.28563 

8999.5604 

- 

- 

7 

174.53673 

1435.5827 

38236.779 

- 

8 

102.91233 

550.32429 

5264.9805 

161461.47 

9 

70.867697 

291.43159 

1792.5214 

19552.206 

10 

53.636458 

183.70287 

860.66613 

5984.7415 

11 

43.221713 

129 .09479 

499.58665 

2626.7294 

12 

36.397896 

97.654760 

327.19911 

1411.8412 

13 

31.656067 

77.869091 

232.87672 

865.00224 

14 

28.210889 

64.574614 

176.05489 

580.69030 

15 

25.619289 

55.182920 

139 .28851 

416.88964 

16 

23.614981 

48.281689 

114.15301 

314.98496 

17 

22.029565 

43.047217 

96 . 203858 

247.69855 

18 

20.751990 

38.972527 

82.926624 

201.11392 

19 

19.706414 

35.731426 

72.816049 

167.59491 

20 

18.839470 

33.106073 

64.927862 

142.69602 

21 

18.112641 

30.946312 

58.645964 

123.69924 > 

22 

17.497498 

29.145784 

53.554692 

108.87271 

23 

16.972648 

27.627345 

49.365398 

97.073311 

24 

16.521712 

26.333870 

45.872543 

87.523269 

25 

16.131961 

25.222272 

42.926505 

79.679145 

26 

15.793363 

24.259519 

40.416223 

73.152296 

27 

15.497913 

23.419916 

38.257813 

67.658992 

28 

15.239157 

22.683208 

36.386929 

62.988249 

,29 

15.011834 

22.033240 

34.753512 

58.980572 v 

30 

14.811617 

21.456984 

33.318124 

55.513588 v 

31 

14.634920 

20.943833 

32.049334 

52.492136 

32 

14.478744 

20.485070 

30.921827 

49.841302 

33 

14.340564 

20.073475 

29.915018 

47.501453 

34 

14.218244 

19.703018 

29.012018 

45.424631 

35 

14.109959 

19.368630 

28.198859 

43.571903 

36 

14.014148 

19.066020 

27.463399 

41.911387 

37 

13.929463 

18.791535 

26.797367 

40.416764 

38 

13.854738 

18.542048 

26.191010 

39.066148 

39 

13.788959 

18.314864 

25.637816 

37.841211 

40 

13.731238 

18.107655 

25.131790 

36.726506 

41 

13.680796 

17.918395 

24.667782 

35.708935 

42 

13.636950 

17.745316 

24.241344 

34.777329 

43 

13.599093 

17.586870 

23.848620 

33.922115 

44 

13.566688 

17.441691 

23.486247 

33.135045 

45 

13.539257 

17.308577 

23.151285 

32.408983 

46 

13.516373 

17.186460 

22 . 841152 

31.737722 

47 

13.497654 

17.074390 

22.553569 

31.115847 

48 

13.482758 

16.971523 

22.286524 

30 .538611 

49 

13.471377 

16.877101 

22.038227 

30.001839 

50 

13.463232 

16.790445 

21.807086 

29 .501845 
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APPENDIX VI 

The Normal Component of the Incident Magnetic 
Field on the Sunlit Side of the Moon 


Our aim is first to express the normal component of 
the incident magnetic field on the sunlit side in a 
Fourier series of the § coordinate and then to truncate 
this series according to a criterion which permits us to 
estimate the resulting error. 

From Equation (2.2.1)/ we can express this component 
as follows: 


H«a^ = [H ( sin0 sincf)costy - cos0sini|/) 


+ H sin0cosd)]e 
ox 


iiJyR^cose + ctsing>) 


where 


a = k. R-.sinQ 
J_ M 


If we use the following identities (see, for ex- 
ample, Morse et al, p. 620) 


e iasin* = £ e im<t» (a) 

nv=-ca 


We can rewrite the above expression as follows: 


H-a 


r 



( Z a cosm4>+ I b sinm^) 
m= 0,2,4.? m=l , 3 , 5 


e 


ik // R M COS0 


ik // R cos6 

+ H ( - cosiwfr + Z, d sinm<*>)e A 

m=l ,‘3.5. . ' ' m~ 2 , 4 , 6 . . m 



where 


a m = " ib m = icos0cos^[J n(+1 (a)-J m _ 1 (a)] 


-e^cosecos^J^ (a) 

c m = " id m = sin9 [ J m-l (a) + J m+l (a)] 
and e m is the Neumann factor [see Equation (2.4.1)]. 


The following rules were used to truncate this series 

1. The terms characterized by m = 1 are always 
kept and if H Q ^ 0, the term m = 0 is also kept. 

2. At a given frequency, we keep all the terms 
in the summation specified by a value of m such that 

«y < . 


x m-l 
(m-1) • 


> 0.01 


where 


X = 


k i R M 


It should be pointed out that in this work the 
parameters of the source were restricted to the following 
values: v sw - 200 km/sec and i p f 60°. Moreover, since 

the frequencies of interest are smaller than 0.04 Hz we ob- 
tain X - 1.9 = y . This value of y insures that the maximum 
value of m that can be kept in the summation is m - 8, 
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and also that if an harmonic specified by m (m > 1) is 
kept then the harmonic specified by m-1 is also kept. 

Our criterion also permits us to establish an upper 
bound of a few percent on the relative mean square error 
resulting from the truncation. In order to show this, 
we first establish two useful inequalities: 

1. Let M be the smallest value of m (9 £ m > 2) 
for which 


(M-1) l 


£ 0.01 


Then 


£ - 0.235 = a 
M 

This result is established by evaluating 
1/{M-1) 

[O.Gl(M-l) I] 7 


M 


for M = 2 to 9. We find that the 


maximum value of this quantity is reached for M = 9 and 
is given by the value of a above. 

2. From an inequality satisfied by cylindrical Bes- 
sel function (see Watson, 1966, p. 16) we have 

J ta) < (g^r -1 a 2 /4m 

J m-1 (0 ° “ (m-1) ! e 

If a = 2Xsin0 and m = M + n (where M and X are defined 
above) , we obtain 

.M-1 


X 


J m-l (a) " (M-1)! 


n MU n m c n 

j e = 0.016 a 


If we assume that H q = 0 and used the above inequalities, 
the relative root mean square error can be bounded as fol- 


lows : 
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sine 


00 

: £ 

m=M 


2 2 
(c m +d m> 1 


1/2 


= [ E ^ ra -l (a > +J m +1 ( “ )}1 
m-M 


1/2 


< 0.016 [ Z 

n=o 


, n n+2 , 2 . 
(a +a ) j 


1 / 2 - 


(1+a 2 ) 

0.016 .. 2,1/2 
(l-o ) 


* 0.017 

So the truncation involves a relative error generally less 

than 2%. The same result can be obtained when H / 0. 

ox 

The above criteria were also used to truncate the repre- 
sentation of the field inside the moon and in the cylindrical void. 

The truncation errors in our calculation were of the same order as those 
derived above. A typical result is illustrated in Figure AP-VI. 

This is a plot of the absolute value of the partial field 

& 

*■ 

associated with each harmonic (m = 1 - 6) of the normal 
component of the poloidal H magnetic field at the lunar 
surface. The parameters of the source are H q = 0 , H qx = 1, 
i p = 60° , Vg W = 300 km/sec and the frequency is 0.04 Hz.' 

In this case, our criteria require that only the harmonic 
m = 1 to 6 be kept in the summation. Actual calculation of 
the term m = 7 reveals that its inclusion would modify 
the. total field by no more than 0.5%. 

For given parameters of the source, the number of 
harmonics kept in our representation will depend on the 
frequency. An example of this is shown in Table AP-VI, 
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where the source used was the same as the one discussed above. 


Table AP-VI 
Fourier Harmonic 


Frequency (Hertz) 
0.0002 
0.0005 
0.005 
0.02 
0.04 


Number of Harmonics 
1 
2 

3 

4 
6 
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Absolute value of the partial field associated with each Fourier 
Harmonic of the radial component of the magnetic field vs. Polar 
Angle, for m - 1, 6. 
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